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Much of our understanding of vibrational excitations and elasticity is based upon analysis of 
frames consisting of sites connected by bonds occupied by central-force springs, the stability of 
which depends on the average number of neighbors per site a. When z < Zc ~ 2d, where d is 
the spatial dimension, frames are unstable with respect to internal deformations. This pedagogical 
review focuses on properties of frames with 2 at or near Zc, which model systems like randomly 
packed spheres near jamming and network glasses. Using an index theorem, Nq — Ns = dN — Nb 
relating the number of sites, N, and number of bonds, Nb, to the number, Nq, of modes of zero 
energy and the number. Ns, of states of self stress, in which springs can be under positive or negative 
tension while forces on sites remain zero, it explores the properties of periodic square, kagome, and 
related lattices for which z = Zc and the relation between states of self stress and zero modes in 
periodic lattices to the surface zero modes of finite free lattices (with free boundary conditions). It 
shows how modifications to the periodic kagome lattice can eliminate all but trivial translational 
zero modes and create topologically distinct classes, analogous to those of topological insulators, 
with protected zero modes at free boundaries and at interfaces between different topological classes. 


I. INTRODUCTION 


Understanding and controlling mechanical stability is important to fields ranging from structural engineering to 
granular materials and glasses. We do not want buildings or bridges to fail, and we want to be able to manage the 
elastic response of materials of everyday life. This review focusses on the elastic and dynamical properties of periodic 
ball-and-spring networks that are at or near mechanical collapse. This may seem like a fairly narrow subject, but it 
is one of considerable richness and impact. 

In a remarkable 1864 paper [1], James Clerk Maxwell undertook the first systematic study of the mechanical stability 
of frames consisting of points, which we will refer to as sites, with connections, which we will usually refer to as bonds, 
between them as a model for such real-world structures as the Warren Truss (patented in 1848) shown in figure 1. 
He defined a “stiff’ frame as one in which “the distance between two points cannot be altered without changing the 
length of one or more connections”. He showed that a stiff frame containing N sites in d dimensions requires 


NB=dN- f{d) 


( 1 ) 


connections, where f{d) = d{d+l)/2 is the number of rigid translations and rotations under free boundary conditions. 
Under periodic boundary conditions, which Maxwell did not consider, f{d) = d. This relation, often referred to as 
Maxwell’s rule, can be reexpressed as a critical coordination number {z = 2 Nb/N), 
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If 2 < z^, the system is not stiff, and if z > z^, it is. As we shall discuss in more detail in the next section. 
Maxwell’s rule for the stability of frames requires modification [2]. Nevertheless, it provides a useful and universally 
used benchmark for the analysis of the stability of frames. We will refer to free frames, i.e., ones under free boundary 
conditions, satisfying Maxwell’s rule as Maxwell frames. It is fairly common practice to use the term isostatic for frames 
satisfying Maxwell’s rule. Though isostatic frames do satisfy Maxwell’s rule, they have more restrictive properties, 
which we will discuss in section H. 

Our principal interest here is in frames whose sites can be collected into identical contiguous unit cells whose origins 
lie on a Bravais lattice and which fill some region of space. We will refer to these frames as lattices and if they are 
subjected to periodic (free) boundary conditions as periodic (free) lattices. Any frame, even those that are not lattices, 
e.g., whose sites are randomly distributed, can also be subjected to periodic (free) boundary conditions in which case 
we will refer to them as periodic (free) frames. For reasons that we will justify more fully later, we will use the term 
periodic Maxwell frame (lattice) for periodic frames (lattices) with average coordination number Zc = 2d = zf° rather 
than z^. Free frames can be liberated from periodic ones by cutting of order bonds. 

Maxwell’s analysis applies to frames with an arbitrary number of sites and bonds. In addition to being a workhorse 
of the structural engineering community [3, 4], it has seen extensive use (though not always with attribution to 
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FIG. 1: (a) The Warren Truss. This is an isostatic structure composed of equilateral triangles with N = 9 sites and Nb = 15 
connections, so that 2N — Nb = 3. The lower left site (indicated by the small triangle) is fixed with respect to the earth, and 
the lower right site is constrained to move only horizontally along a track with wheels. These constraints reduce the number 
of free degrees of freedom of the sites by 3 to Wree = 2N — 3 = 15 = Nb- (b) A reduced version of the isostatic Warren Truss 
indicating conventions for labeling sites and bonds. The arrows indicate indicate bond directions following the convention 
described in section IIB, Eqs. (12) to (15). 


Maxwell) in physics, materials science, and mathematics. It is a critical component of the theory of structural glasses 
[5-9], rigidity percolation [10-14], framework silicates like /3-cristobalite [15], jamming of packed spheres [16-20], 
biopolymer networks [21-28], and of some theories of protein folding [29]. Rigidity percolation and jamming generally 
involve central forces only, and the Maxwell relation and its generalization can be applied to these problems directly. 
Structural glasses and biopolymer networks have bending forces that require a modification of the Maxwell rules, and 
we will not make much contact with them in what follows. 

Both randomly packed spheres (figure 2(a)) and diluted elastic networks (figure 2(b)) pass from a state that does 
not support stress to one that does as the number of bonds or contacts increases. At the critical point separating these 
two states, the coordination number is at or near the Maxwell critical value of Zc- The rigidity percolation transition 
is generally continuous [10-12, 14] in two-dimensions, and it is well-described by the language of critical phenomena, 
applied so successfully to percolation in random resistor networks [30, 31]. Both shear and bulk elastic moduli increase 
from zero as a power law in (z — Zc), and there is a divergent length scale associated with the probability that two 
sites are in the same rigid cluster as a function of their separation. In three dimensions, the rigidity transition is 
apparently first order [14]. 

Jammed systems exhibit different behavior. They are constructed by increasing the density of spheres in a fixed 
volume until they first resist a further increase. At this point, they are jammed, and their bulk modulus B (which 
resists compression) is nonzero [32, 33], but their shear modulus G is of order 1/N [34, 35]. At large N, G increases 
linearly in Az = (z — Zc) [32, 33, 36]. The density of states of systems with Az > 0 exhibits a crossover from 
Debye-like behavior (^ in d dimensions as a function of frequency w) at low frequency to a flat-plateau beyond 
a characteristic frequency uj* ^ (Az) [33, 37]. There are two diverging length scales, I* ^ (Az)“^ and It ~ (Az)“^/^, 
which can be extracted [37] by, respectively, comparing uj* to the longitudinal and transverse sound frequencies clI~^ 
and ctI~^ ~ (Az)^/^Z“^ at wavenumber q ^ 1“^, where cl and ct are the longitudinal and transverse 

sound velocities. Other interpretations of these lengths invoke finite-size effects in the modes of a finite sample cut 
from a larger one [17, 18, 38] for I* and stability of the system to boundary perturbations [39, 40], crossover between 
strong and weak scattering at w* [41, 42], or correlations in systems with z < Zc [43, 44] for It- 

Though the transition to elastic stability in both rigidity percolation and jamming occurs at or near the Maxwell 
critical point z = Zc, properties at and above the two transitions are quite different. Presumably, the difference 
is a reflection of the different geometrical arrangements of the Maxwell lattices at or near the two transitions. An 
interesting question then is what precisely are the differences. Can they be quantified, and if so how? These questions 
then raise a broader question of whether there are wider classes of Maxwell frames that lead to elastic and vibrational 
structure different from those of the percolation and jamming lattices. As a first step toward answering that question, 
it is useful to study periodic Maxwell lattices, like the square and kagome lattices shown in figure 3 with NN bonds 
and z = 2d = 4, and the reduced coordination lattices they spawn when bonds are cut to create free surfaces. They 
can easily be moved away from the Maxwell limit toward greater rigidity by uniformly [45] or randomly [46] adding 
further neighbor bonds, and they have an advantage that they lend themselves to exact calculations of elastic response 
both in the uniform and random systems (via effective medium theory). In addition, in spite of their simplicity, they 
present a surprisingly rich phenomenology that inform us generally about Maxwell frames. 

Much of the language [2, 47-49] for probing the mechanical properties of frames used in this review was developed by 
members of the Structural Mechanics group at the University of Cambridge, which provided an elegant generalization 
of Maxwell’s relation based on general principles of linear algebra and used it to deepen our understanding of the 
more subtle properties of frames. Though well-known in the engineering community, it is less so in the physics and 
materials science communities. This review will, therefore, begin in section II with a fairly comprehensive review 
of this work. It is followed by five more sections and four appendices. Section III deals the long-wavelength elastic 
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(a) (b) 


FIG. 2: (a) packed discs of two sizes just above the jamming transition. The dark red lines are chains of force that are a 
response to the pressure required to pack the particles at ^ > 2 c (courtesy of Carl Goodrich), (b) A representative bond-diluted 
lattice near the rigidity percolation threshold of ~ 3.96. 



(a) 



FIG. 3: (a) Square, (b) distorted square, (c) kagome lattice, and (d) twisted kagome lattice 


limit and how it can be calculated using the developments in references [47, 48]. Section IV tailors the formalism 
of the preceding sections specifically to periodic lattices. Section V looks in detail at the properties of the three 
simple lattices in figure 3, and section VI explores lattices with topologically protected states at interfaces, which are 
closely related to topologically protected electronic states in the quantum Hall effect [50, 51], polyacetylene [52] and 
topological insulators [53-58]. Section VII presents some final thoughts and speculates about future directions. The 
four appendices provide mathematical detail and display derivations of various important relations. 

Section V contains several subsections. Section V A with its simple analytical treatment of the square lattice sets 
the stage for the study in section VB of the simple kagome lattice. Both lattices have lines of zero modes in their 
Brillouin zones arising from their having straight sample-traversing filaments of colinear bonds. Section V C then 
explores how the simple geometrical operation of “twisting” neighboring triangles in the standard kagome lattice 
]figure 3(c)] to convert it to the twisted lattice [figure 3(d)] gaps the phonon spectrum of the simple kagome lattice 
at all wavenumbers except at the origin and leads to zero-frequency Rayleigh surface modes [59] not present in either 
the square or standard kagome lattices. 

This review is intended to be as much a pedagogical introduction to periodic Maxwell lattices as an overview of 
the subject. Except in section VI, which requires the use of some fairly subtle concepts detailed in reference [60], 
it provides sufficient calculational detail that even someone totally new to the subject should be able to follow it. 
Though, as the introductory paragraphs have indicated, this review was inspired in part by jamming and rigidity 
percolation, it is not about these subjects, and they will be considered only when they have direct overlap with the 
ideas being presented. 
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FIG. 4: (a) to (c) Frames satisfying the Maxwell rule, (a) has 6 sites, 7 bonds, 5 zero modes, and two mechanisms indicated by 
the dotted bonds, (b) has 6 sites, 8 bonds, 4 zero modes, and one mechanism, (c) and (d) are constructed from (b) by adding 
an additional diagonal bond, (c) satisfies the Maxwell rule with only the three trivial zero modes, (d) has 4 zero modes and 
one state of self stress indicated by the arrows on the bonds in the left square. 


II. GENERALIZED MAXWELL RELATION AS AN INDEX THEOREM 
A. The Maxwell rule and states of self stress 

Each site in d-dimensions has d independent translational degrees of freedom, and in the absence of constraints on 
point motion, a collection of N points without connections has Nf^ee = dN zero-frequency displacement modes, which 
we will refer to as zero modes. In the presence of constraints, will be less that dN. Each connection reduces the 
number of zero modes by one. Thus if there are Nb connections and no constraints, there are 

No = dN- Nb (3) 

zero modes. Of these, /(d) are the trivial ones associated with rigid translations and rotations. Any other zero modes 
involve internal displacements of the sites and are generally called mechanisms [2] in the engineering literature and 
floppy modes in the physics literature [6]. Equation (3) reexpressed in terms of the number of mechanisms M is 

M = dN - Nb - fid). (4) 

We will refer to Eqs. (3) and (4) as the Maxwell count. A frame is stiff if it has no mechanisms. Setting M = 0 yields 
the Maxwell rule [Eqs. (1)]. Eigure 4(a) depicts a simple frame that obeys Maxwell’s count. It consists N = 6 sites 
and Nb = 7 bonds, and it has Ag = 2x6 — 7 = 5 zero modes and M = Nb = 5 — 3 = 2 mechanisms. 

The simple Maxwell rule does not apply to all frames [2]. Consider the two-square frame with N = 6 sites and 
Nb = 8 bonds shown in figure 4(b). It has one mechanism as expected from the Maxwell count. If an extra bond is 
added. Maxwell’s rule would say that the frame is stiff with no mechanisms. The extra bond, however, can be placed 
as a diagonal in the right square [figure 4(c)] or as an extra diagonal in the left square [figure 4(d)]. In the first case, 
there are no mechanisms, and Maxwell’s rule applies. In the second case, however, the mechanism present before the 
extra bond was added remains, and the Maxwell count is violated. But the left square with crossed diagonal bonds 
has an extra redundant bond not needed for its rigidity. It also has a new and interesting property: the outer bonds 
of the square can be placed under tension (compression) and the inner diagonal bonds under compression (tension) 
such that the net force on all sites is zero. This is a state of self stress, which, because of its repeated use in this 
review, we will usually abbreviate as SSS. This theme can clearly be repeated with each added bond either decreasing 
the number of zero modes or increasing the number of states of self stress to yield the modified Maxwell count [2]: 

No = dN-NB + Ns or No-Ns = dN-NB, (5) 

where Ns is the number of SSSs. This is an index theorem [61, 62] , which we will derive in section IIB, relating 
mode and self-stress count to geometric properties of the lattice. We will refer to it simply as the Index theorem . 

Two types of mechanisms can be distinguished: “finite” ones in which finite-amplitude displacements of sites stretch 
no bonds and “infinitesimal” ones in which bond lengths do not change to first order in the magnitude of displacements 
but do so to second (or higher) order. The Index theorem [2], as we shall show below, follows from the assumption of 
a linear relation between site displacements and bond lengths, and it applies only to infinitesimal displacements, i.e., 
it counts both finite and infinitesimal mechanisms but does not identify which is which. Eigures 5(a)-(b) show how a 
finite mechanism can be converted into two infinitesimal mechanisms and one SSS. A configuration of self-stress that 
is particularly important for the current study is any straight line of bonds under periodic boundary conditions, which 
we will refer to as straight filaments, as shown in figures 5(c)-(d). Changing the straight filament to a zigzagged one 
removes this state of self stress. On the other hand the “zigzaging” periodic ladder configuration shown in figure 5(e) 
has one SSS, rather than the two that a straight ladder would have. Tensions alternate in sign from bond to bond 
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FIG. 5; States of self stress: (a) A frame with four sites (1-4) and four bonds (italic 1,2,3, and 4) with 4 zero modes and one 
finite mechanism (dashed lines); (b) A frame in which the length of bond 4 is equal to the sum of the lengths of bonds 1 to 3. 
Now there is a SSS in which bond 4 is under compression (tension) and the other three are under tension (compression). Both 
sites 2 and 3 can undergo infinitesimal displacements without changing the length of any bonds, and there are two infinitesimal 
mechanisms, (c) A line of parallel bonds forming a sample-traversing filament under periodic boundary conditions as depicted 
in (d). The forces on all sites are zero if all of the bonds are under equal tension or compression, and there is one SSS . (e) A 
zigzag ladder under periodic boundary conditions with one SSS . 


in this SSS, a property, which will be important in what follows, that prevents it from having any zero wave-number 
component. 

A system in which there are neither any mechanisms (M = 0) nor any states of self-stress {Ns = 0) is isostatic . 
A finite isostatic system necessarily satisfies the Maxwell relation z = , but a system with z = z^ can have any 

number of mechanisms provided it is equal to the number of SSSs. The distinction between satisfying the Maxwell 
rule and being isostatic is often lost in the literature, and it is common practice to refer to any system that satisfies 
Maxwell’s rule as isostatic. In this review, we will keep the distinction, referring to any free frame satisfying Maxwell’s 
rule as a Maxwell frame, reserving the term isostatic for those free Maxwell frames satisfying M = Ns = 0. As we 
shall see in section IIIC, the extension of this definition to periodic frames presents some problems [49]. Since the 
term isostatic has become so prevalent, we propose in that section a definition of this term that is in the spirit of the 
definition for free frames and consistent with common usage for periodic frames. 


B. Equilibrium and compatibility matrices 

In the absence of external forces, the equilibrium force at each site in a frame is determined by the tensions in the 
bonds it shares with other sites. This is true whether or not the site is in mechanical equilibrium; if the force at a 
site is nonzero, the mass at that site will accelerate according to Newton’s laws. If forces at sites arising from bond 
tensions are nonzero, they can be balanced by external loads to create an equilibrium situation in which the total 
force on each site is zero. Clearly, in mechanical equilibrium, the external loads are the negative of the forces at each 
site arising from bond tensions. 

For central forces, the tension in a bond is parallel to the bond. Thus its direction is specified by bond orientation, 
but its magnitude and sign can vary. Let F be a vector in the diV-dimensional space, Vp) of the d components of 
force at each site on the lattice exerted by tensions in the bonds that terminate on it, and let T be a vector in the 
iVs-dimensional space, Vj, of the of the bond tensions, which are scalars of either sign. External loads at sites are 
represented by the diV-dimensional vector L. In equilibrium when sites do not accelerate (or move if there is external 
friction), L = — F. Since the relation between forces and tensions is linear, there is a dN x Nb dimensional matrix Q 
(with dN rows and Nb columns), called the equilibrium matrix, that maps Vj to Vp: 

QT = -F = L, (6) 

where the hnal relation only applies in static situations. 

The null space or kernel of Q, ker(Q) of dimension nullity(Q), is the set of all vectors mapped to the zero vector. 
Any vector in the null space of Q represents a state of self stress because it corresponds to tensions on a set of 
bonds for which the forces on all sites are zero. Thus nullity(Q) is equal to the number of SSSs Ns- Vectors T not 
mapped into the null space of Q are in the orthogonal complement. Opt. Commun. (Q), of ker(Q). The dimension of 
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Opt. Commun. (Q) is equal the to the rank of Q. The rank-nullity theorem of linear algebra [63] relates the rank and 
nullity of a matrix to its column dimension: 

rank(Q) -|- nullity(Q) = rank(Q) -|- Ns = Nb- (7) 

Elongation of bonds are determined by the displacements of the sites to which they are attached. The elongations of 
individual bonds are necessarily parallel to bond vectors for central forces. Let U be a vector in the dfV-dimensional 
space, Vu, of site displacements and E be a vector in the TVs-dimensional space, Ve, of bond elongations. The 
iVs X c?A^-dimensional compatibility matrix C maps Vy to Ve: 

CU = E. (8) 

The null space of C is the set of displacements U that do not change the length of bonds, i.e., the set of zero modes 
of the system; thus, nullity(C) = Nq. The rank-nullity theorem applied to C yields 

rank(C) -|- nullity(C) = rank(C) -|- A^b = dN. (9) 

The equilibrium and compatibility matrices are not independent: they are matrix transposes of each other. To see 
this, we can calculate the work done under infinitesimal distortions of the system in the presence of sites forces (and 
thus necessarily tensions in the bonds) in two ways: First the work Wl done by external loads in displacing sites and 
second the work Wt done by bond tensions in stretching bonds. The two calculations must yield the same result: 

Wl = L^U = T^Q'^U = Wt = T'^E = T^CU, (10) 

where the superscript T refers to the transpose of a matrix. Since this relation is valid for all U (even U in the null 
space of C) it must be that C = Q^. The rank of a matrix is equal to the rank of its transpose, rank(Q) = rank(C), 
and subtracting (7) from (9) yields the Index theorem of (5). The equilibrium and compatibility matrices have proven 
useful in many contexts, and we mention only one in which frames with spring lengths that do not correspond to the 
length of the bonds they occupy are constructed from frames in which they do [64]. 

To construct Q and C for a particular frame, it is necessary to have explicit representations for the vectors U, E, T, 
and F. A frame consists of N sites labeled s = 1, ...,N at equilibrium positions R{s) connected by Nb bonds labeled 
/3 = 1,..., Nb- Under frame distortions, site positions change to 

X{s) = Ris)+u{s), (11) 

where u{s) with components Ui{s), i = is the displacement vector at site s. The force at site s is f{s) = 

(/i(s), .../d(s)), and the dA^-dimensional displacement and force vectors are, respectively, U = {ui, ...,un) and F = 
(/i, /iv)- 

Each bond /? = [s,g, s^] connects a pair of sites sp and s]^, whose separation in the equilibrium frame is the vector 

HP) = Hlsp, 4]) = -^(4) “ = 44. (12) 

from Sjs to s'^, where is the length of bond /3 and bp is the unit vector along bond j3. The arrows in figure 1(b) 

show an arbitrarily chosen choice of directions of vectors bp for bonds in a simple frame. Other choices will lead to 
different equations relating forces to tensions, but no the physical tensions in the frame. Let the tension in bond (3 
be tp, and associate with that bond a vector tension, 

tp = tpbp. (13) 

With this convention, the force exerted on sites sp and s'p by bond [s/3, 4] 9-re, respectively, tp and —tp, or, equivalently, 

the force on a site s from a bond /3 that it shares with another site is tp if bp points away from the site and —tp 
if it points toward the site. With these definitions, we can construct the compatibility matrix from the the bond 
elongation relations, 

ep = bp- {u{s'p) - u{sp)), (14) 

and the equilibrium matrix from the site force equations, 

Zs 

fH) = £ sign(/3)t/3, 

/3=i 


( 15 ) 
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where Zg is the number of bonds site s shares with its neighbors and sign(/3) = +1(—1) if the arrow of bond /3 points 
away from (toward) site s. 

It is instructive to carry out explicit calculations of Q and C for a simple frame. We consider the frame shown in 
figure 1(b) with = 4, Nb = 5 and 2N — Nb = 3, the number of modes of rigid translation and rotation. This 
frame has no floppy modes and no state of self stress, and it is isostatic. The hgure labels sites (s = 0, ...,3) , bonds 
(/3 = 1,..., 5), and bond directions. The five displacement-elongation equations are 

ei = 6i • (mi - Mo) 62 = ^2 • {U2 “ Mi) 63 = 63 • (m2 - M3) 

64 = 64 • (m 3 - Mo) 65 = bs • (m 3 - Ml), ( 16 ) 

and the four vector force equations are 

/o = + ti fi = —ti +^2+^5 

/2 = —^2 — *3 /a = ^3 — ^4 — ^ 5 - ( 17 ) 

The 8x5 compatibility and 5x8 equilibrium matrices are easily constructed from equations ( 16 ) and (17). If, however, 
we are interested only in internal deformations of the frame and not its uniform translations and rotations, we can 
apply three constraints to the motion, most easily by pinning site 0 so that Mo = 0 and placing site 1 on a horizontal 
rail as shown in figure 1(a) to fix ui^y = 0 so that A^free = 5. We also allow /o and fi^x to take on whatever values 
needed to satisfy the constraints, so they do not enter into our equations. This leaves us with a 5 x 5 compatibility 
matrix, 


C = 


/ 1 0 0 
' _1 1 ^ 
2 2 2 
0 1 0 

0 0 0 

V1 0 0 


0 

0 

-1 

1 

2 

_ 1 
2 


0 

0 


(18) 


mapping U = {ui^x,U 2 ^x)U 2 y,u^^xi'u^,y) to E = (61,62,63,64,65). The equilibrium matrix, mapping T = 
(ti,t 2 , 13 , 14 , 15 ) to L = -F = -{fi,y, f 2 ,x, f 2 ,y, f 3 ,x, f 3 ,v) Constructed from Eq. (17), is trivially equal to C"^. Both 
Q and C are square invertible matrices: their nullspaces are empty, both A^o = 0 an'l M = 0, and the system is 
isostatic as required. Thus, the tensions on the bonds are uniquely determined by the forces on the sites and vice 
versa: the frame is statically determinate. And the elongations of the bonds are uniquely determined by the site 
displacements and vice versa: the frame is kinematically determinate. Thus, an alternative, and perhaps preferable, 
definition of an isostatic frame is that it be both statically and kinematically determinate. Another way of dealing 
with the trivial zero modes is to introduce “reaction forces” to yield 8x8 matrices Q and C. 


C. The dynamical matrix 

So far, we have only discussed tensions and stretches of bonds without specifying any relation among them. In our 
“ball-and-spring” frames, each bond (3 is occupied by a Hooke’s law spring whose energy is half its spring constant kb 
times the square of its elongation. Let k be the Nb x Nb diagonal matrix of spring constants, then the elastic energy 
of the lattice is 

Eei=iE^kE = ^U^KU, (19) 

where 

K = QkQ^ = C'^kC (20) 

is the dN x dN stiffness matrix. Normal-mode frequencies depend on mass as well as the stiffness matrix. The kinetic 
energy requires the introduction of a mass matrix M. We will restrict our attention of frames in which the mass of 
all mass points is equal to m, in which case M = m I, where I is the unit matrix, and the kinetic energy is 

Akin = imU^U, (21) 

where U is the velocity vector. Normal modes are then eigenvectors of the dynamical matrix: 


( 22 ) 



The Lagrangian in the presence of eternal loads is thus 


L = imU^U - Fei - U^L, 


(23) 


and the equation of motion is 


mU = -^-L = -KU-L = F-L, ( 24 ) 

which vanishes when the external load L is equal to the force F = — KU = —QT exerted by bond stretching. Note 
that the equilibrium matrix can be used to calculate F whether or not the system is in static equilibrium or not. On 
the other hand L = QT only in equilibrium when there is no acceleration (assuming no friction forces). 


III. THE ELASTIC LIMIT 


A. Strain and the elastic energy 


Strain is a measure of macroscopic distortions of an elastic medium. The macroscopic deformation tensor A = J + jy, 
where I is the unit tensor, determines displacements at boundary sites ss of either a finite sample or the periodic 
box under periodic boundary conditions: 

X{sb) = >^R{sb), ot u{sb) = vRisB), (25) 


and the macroscopic strain tensor is 

£ = i(A^A-/) = + + « U'^ + V^), ( 26 ) 

where the final form is the linearized limit, which is all that concerns us here. 

The elastic energy density associated with the macroscopic strain is 

fel 

where Kijki is the elastic constant tensor and the Einstein convention on repeated indices is understood. The elastic 
strain Sij is symmetric and has Ud = d{d+ l)/2 independent components in d dimensions. It can be expressed [65, 66] 
as an Od-dimensional vector (Voigt notation), which in two dimensions takes the form. 


— \^xx^ ^yy^xy) ’ 

The elastic tensor is then an ad x ad matrix, which in two dimension is 


/ ^xxxx 

K 

^^xxyy 

‘2K 

^^^xxxy 

IK = 1 Kxxyy 

^yyyy 

2K 

^-^^yyxy 

\^^xxxy 

2K 

^-^^yyxy 

4K 

^^^xyxy 


isotropic 


[b + g b-g o\ 

\b-G B + G q \ , 

\ 0 0 Gy/ 


(27) 


(28) 


where the final form is the isotropic limit, where B = A + /i is the bulk modulus and G = /i is the shear modulus with 
A and /i the standard Lame coefficients [59]. The linearized Cauchy stress tensor is 

^ij — XijkGkl- (29) 


Mechanical stability requires that all ad eigenvalues of the Voigt elastic matrix K be positive. Thus the elastic 
energy of an elastically stable system can be expressed as a sum of squares of ad independent linear combinations of 
strains (the eigenvectors of the elastic matrix) with positive coefficients (the eigenvalues of the elastic matrix). We 
shall see shortly that some or all of the eigenvalues of the elastic matrix in lattices at or near the Maxwell limit may 
be zero. 


B. Elastic limit and States of Self Stress 

Calculations of the elastic tensor require some method of maintaining strain. The usual picture is that boundary 
sites of a finite frame are clamped to conform with (25). Since these sites are fixed, their displacements and the forces 
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associated with them do not enter the calculation of Q and C. An alternative approach, which we implement, is to 
apply periodic boundary conditions (PBCs) to the frame, which may or may not be composed of repeated unit cells. 
In this approach, it is the boundaries of the periodic cell that satisfy (25). 

If the positions of all sites, rather than just boundary sites, displace according to Eq. (25), the elastic distortion is 
affine. Under such distortions, the relative displacement of the sites associated with bond /? is 

u{sp) - u(s^) = T] - Risp)) =ribp, (30) 


and the affine stretch of bond (3 is 


ef = bl\bp= (31) 

where we used the fact that bpffipj is symmetric in i and j to convert ?7 to e in the linearized limit. The affine elastic 
energy density is then 


/ef = (32) 

where V is the volume and EaS is the vector of affine elongations 

Affine response throughout a sample is the exception rather than the rule. It is guaranteed to occur in absolutely 
homogeneous systems and in lattices with one site per unit cell, but it can occur in certain other systems with special 
relations among elastic constants or special arrangements of sites in periodic systems with multi-site unit cells (for 
example, as we shall see, in the kagome lattice). Generally, however, the forces at at least some sites under an affine 
strain imposed by macroscopic stain at the boundary are nonzero, and these sites will relax to positions of zero force 
(or equivalently until the energy reaches its minimum value), in which cases their displacements are nonaffine. In A, 
we show that the result of this relaxation is that the elastic energy density becomes [47, 48, 67] 

/el = ^Erff,J(k-i),,)]-iEaff,, .t„)2 (33) 

a 

where Eaff,s and (k“^)ss are the projections of EaS and k~^ onto ker(Q), and ta is the ath orthonormal basis vector of 
ker(Q). Thus, only the projections of the affine displacement vectors onto states of self-stress contribute to the elastic 
energy. 

Equation (33) encodes a great deal of information. 

• First, it shows that lattices cannot be elastically stable unless they have SSSs in the presence of conditions that 
constrain the macroscopic strain - a simple reflection of the fact that forces on each site must be zero once 
equilibrium is reached in the presence of imposed strains, which necessarily induce bond tension. 

• Second, only those states of self-stress with a nonzero overlap with the affine bond elongations contribute to the 
elastic energy. These states necessarily traverse the sample, and they are load-bearing. The straight filament of 
figure 5(c) (wound to a circle - figure 5(d)), whose bonds all have the same sign of tension provides an example 
of a load-bearing state, whereas the zigzag state of figure 5(e) and the localized crossed square of figure 4(d) 
does not. 

• Third, because it is a sum of squares of linear combinations of strain, it shows that there must be at least 
Od = d{d l)/2 load-bearing SSSs to produce an elastically stable system with an elastic matrix with 
positive eigenvalues. 

C. Isostaticity and periodic boundary conditions 

In section II A, isostatic lattices were defined as ones that are both kinematically (Ag = f{d)) and statically 
determinate {Ns = 0). This definition is unambiguous for finite free frames. It would seem natural to define an 
isostatic frame under PBCs in the same way, but there is a problem with this definition [49]. Under PBCs, the shape 
and size of the frame boundary is fixed; the compatibility matrix (and by extension, the equilibrium matrix) applies 
to displacements and does not apply to changes in the shape or size of the periodic boundary [34, 35, 68], which 
are described by macroscopic strain. In order for a lattice to be elastically stable, it must have at least d{d -|- l)/2 
SSSs. Reference [49] defines a lattice under PBCs to be statically determinate if Ns = d{d I)/2 and kinematically 
determinate if Nq = d, but it does not propose applying the terms isostatic to such systems. We propose calling 
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a frame under PBCs isostatic \i Nq = d and Ns lies between 1 and d{d + l)/2. This corresponds more or less to 
common usage in the jamming literature. If greater precision is required, the value of Ns can be specified via the 
term isostatic. 

The above discussion applies to any frame, whether it is a lattice with unit cells on a Bravais lattice or not, and in 
particular to finite-frame approximations to infinite random systems such as those encountered studies of jamming. 
These frames can be subjected to PBCs by wrapping them on a torus to create a single cell in which sites separated 
by repeat distances of the torus are identified. We will refer to these as toroidal PBCs. Since there is only one 
cell with randomly placed sites, there is no wavenumber q to index the dN vibrational states. Alternatively, the 
single (large) cell can be periodically repeated in a Bravais lattice, in which case there are dN vibrational bands with 
wave-number-dependent frequencies. The q = 0 limit of these bands are the vibrational states under the toroidal 
PBCs and are generally the ones of physical interest, though interesting information about the stability of jammed 
structures can be obtained from an examination of spectrum at q ^ 0 [40]. 

Lattices under periodic boundary conditions are special. They consist of N^ periodically repeated unit cells with n 
sites and rib bonds, and the Index theorem reads Nq — Ns = {dn — nb)Nc ■ In a periodic Maxwell lattice, rib = dn, 
and Aj) = Ns- Thus if Aj) has its minimum value of d, Ns = d, and under the definition proposed above, such a 
lattice would be isostatic. It is impossible to have Nq = d and for Ns to have any other value than d in the interval 
between I and d{d -I- l)/2 because a change in Ns requires a change in n or rib and thus of dn — Ub, which would lead 
to a change in Nq — Ns of order Nc^ 1 rather than of order one. Thus we can uniquely define an isostatic periodic 
lattice under PBCs to be one with Nq = Ns = d. 

Finite frames can be constructed from ones subjected to PBCs by cutting bonds along the periodic boundaries and 
“liberating” the sites attached to them. Since opposite sides of the boundary are equivalent, it is only necessary to 
cut bonds along half of the boundary to liberate a full free lattice. Thus, an N^ x Ny free square lattice is obtained 
by cutting N^ -f Ny bonds in the lattice under PBCs. The cutting process thus reduces the number of bonds by of 
order in a lattice of N sites in d dimensions and increases Aq — Ns by the same amount. If the periodic 

lattice is isostatic, then there are necessarily of order extra zero modes in the finite lattice. If on the other 

hand, the lattice under PBCs has of order SSSs that are removed on cutting, there may be no increase in 

zero modes at all upon cutting. 


IV. PERIODIC LATTICES 

Our primary interest is in periodic lattices, and we review here notation that we will use to describe them. 


A. Notation 

A general lattice has Aj, unit cells consisting of n sites and Ub bonds so that the total number of sites and bonds 
are, respectively, N = Nsn an Nb = NcUb- The Bravais lattice vectors Ri, where £ = {li,-- - Id) with each li an 
integer, are linear combinations of the primitive translation vectors aj, j = 1, ■■■ d: 

d 

Ri = ''^ljaj. (34) 

j 

The positions of sites and bonds in cell £ in the undistorted lattice with a basis are 

Re^fj. = Ri + ry,, p = l,---,n 

Re,i3 = Ri + rp, /3 = 1, • • • , nt,, (35) 

where Vy and r,g are, respectively the positions of the sites and bonds relative to the origin of the unit cell. The 
positions of lattice sites in a distorted lattice are 

where Uy{£) with Cartesian components Uyy{£) = u„{£) is the displacement vector of site (£,p)- The components 
of U are thus the dN displacements Ua{£)- The Cartesian components fa{£) of the force vector fy{£) are the dN 
components of F. The Nb bond elongations e^d) and bonds tensions t^d) are the components of the of E, and T, 
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respectively. Fourier transforms in periodic lattices are defined in the usual way in terms of wavenumbers q in the 
first Brillouin zone: 

“a* W = E Mq) = E W (37) 

° q t 

E tp{q) = Y, (38) 

^ q I 

and similarly for f^, ep and other site and bond variables. These quantities can also be defined without the site and 
bond basis vectors, and rp in the exponentials, and we will usually leave them out. They will, however, be of use 
in our discussion of topological states in section VI. The components of the equilibrium matrix Q in periodic lattices 
are Qcr,/ 3 (£, £'). Its Fourier transform is 

Qa/3(g) = E 0). (39) 

i 

Again, the basis vectors in the exponents are not required. The compatibility matrix is C(q) = Q^(<7). 

With these definitions, we can reexpress equations (6) and (8) as 

Q(g)t(9) = -f(g) C(q)u(q) = e(q), (40) 

where the “math boldface” vectors are defined as t(q) = {ti{q), ...,tn^{q)) and u{q) = {ui{q), ...,Un{q)) and similarly 
for e{q) and f(q), and Q(q) is the dn x nt matrix with components Qcrp{q)- There is one equation for each of the 
values of q, giving us as many independent equations as (6) and (8). The Index theorem applies to these equations 
for each q: 

no{q) - n^iq) = dn - rib, (41) 

where no{q) = dimker(C) is the number of zero modes and ns{q) = dimker(Q) is the number of states of self stress 
at wavevector q. Of course, Nq = ’^o(q) and Ns = J2q ^s{q)- 

In periodic Maxwell lattices, dn = nb, and there are exactly as many zero modes as states of self stress at each 
q. Under PBCs, there are always d zero modes that arise from translational invariance, and these are necessarily at 
q = 0, implying that there are at least d SSSs at q = 0. There may be more, but each additional SSS will require an 
additional zero mode, which is a q = 0 mechanism. At nonzero q, there is no general reason for zero modes to exist, 
but if they do, the are necessarily accompanied by SSSs. We will see that this is a common theme in our study of 
specific lattices in section V: removing states of self stress eliminates zero modes and “gaps” the phonon spectrum. 
Following (19), the potential energy in terms of a periodic harmonic lattice is 

Ksi = ^^e^(q)ke(q) = ^ ^ u^(q)K(q)u(q), (42) 

q q 

where k is the nb x nb diagonal matrix of spring constants, and 

K(q) = Q(q) k Q'f (q) = TOD(q) (43) 

is the dynamical matrix. In periodic systems, nearest-neighbor (NN), next-nearest neighbor (NNN), and further- 
neighbor bonds are well defined, and bond vectors can be expressed as the direct sum of NN and NNN components, 
e.g. e = esfN © ^nnn, and elastic constant matrix and dynamical matrices can be decomposed into NN and NNN: 

D(q) = DArAr(q) -|- T)NNN{q). (44) 


B. The elastic limit 

1. The elastic energy 

Under affine strain, the strain of equivalent bonds in different unit cells are identical, and we can describe affine 
strain in terms of the nb dimensional vector eaff,s- As a result, the elastic energy depends only on the projection of 
the affine strain onto the q = 0 SSSs. Following (33) and A, the elastic energy density is thus 

/el = ^e^ff.s(k;/)”^eaff,s ^ 


(45) 
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where Vc = Nc/V is the volume of a unit cell and k^s is the projection of the k on to the nullspace of Q (<7 = 0). The 
final form with {a = 1, ...,ns(0)) basis vectors for the null space of Q (<7 = 0) applies when there is a single spring 
constant k. 

Equation (45) constrains the elastic energy of Maxwell lattices, which depend on the number, mo, of <7 = 0 
mechanisms and on the overlap of q = 0 SSSs with the q = 0 affine bond elongation Gaff. Consider periodic Maxwell 
lattices, for which with Nq = Ns- 

1. mo = 0: In this case, there are exactly d zero modes and exactly d SSSs. There are now two possibilities: 

(a) All d SSSs have a nonzero overlap with the affine bond elongations, and the elastic matrix (28) then 
has d positive and d{d + l)/2 — d = d{d — l )/2 zero eigenvalues, which correspond to zero-energy elastic 
deformations, now referred to as Guest modes [49]. This case corresponds to what we call an isostatic 
periodic lattice. As we shall see, this is the situation for the square, twisted kagome, and topological 
kagome lattices [section VIB] in which there are two positive eigenvalues and one zero-energy elastic mode. 

(b) Fewer than d SSSs have a nonzero overlap with affine bond elongations, and as a result there are fewer than 
d positive eigenvalues to the elastic matrix and more than d{d — I)/2 zero-energy elastic distortions. The 
zigzagged square lattice to be discussed in section VD provides an example of this behavior. Finite periodic 
approximates to infinite-unit cell systems, such as packed spheres at the jamming transition [34, 35] and 
randomized quasi-crystalline Penrose tilings [section VD] [69] exhibit smaller and smaller overlap as the 
order of the approximate increases, leading to shear moduli that vanish as I/n as n —)■ 00 . 

2. mo > 0, and there are d + rriQ SSSs, that may or may not have an overlap with affine bond stretches. If more 
than d overlap, the additional SSSs stabilize elastic response relative that of the lattice with mo = 0, and we are 
presented with the curious situation in which additional zero modes increase elastic rigidity. The usual cause of 
this effect is the appearance of sample traversing straight filaments that support macroscopic stress, but they 
also, as we have seen, introduce additional infinitesimal zero modes. The untwisted kagome lattice, with three 
q = 0 states of self stress produced by parallel filaments along three independent directions, is elastically stable 
with all eigenvalues of the elastic matrix positive. In spite of extra mechanisms it is possible for only some or 
even none of the SSSs to overlap with affine bond elongations, in which case the elastic energy can even fall to 
zero. The latter situation occurs in unrandomized periodic approximates to Penrose tilings [69] in which there 
are of order zero modes and states of self-stress but in which all elastic moduli vanish [section VD]. 


2. Stiffness matrix of 2d periodic lattices with Guest modes 

Two-dimensional periodic lattices with one or two q = 0 SSSs and two q = 0 zero modes have two and one Guest 
modes, respectively. Fully gapped isostatic lattices have two such SSSs , but so do lattices that have additional zero 
modes at non-zero q. When there is only one Guest mode, the elastic matrix of (28), which acts on the strain vector 
ev, has one zero eigenvalue with normalized eigenvector Vq = (uo,i, uo, 2 , "Cos) (be., the strain of the Guest mode with 
amplitude Ug is UgVq) and two positive eigenvalues, Ki and K 2 , with respective associated eigenvectors Vi, and V 2 
so that 


Kij = KiVi^iVij + K2V2,iV2,j (46) 

The long-wavelength stiffness matrix K is determined by K, and as B shows, its determinant depends only on i^i, 
K 2 , and the Guest-mode eigenvector Vg: 

det K(q) = ^ATiATa {vo 2 ql - ‘^vosqxQy + voiqlf- (47) 

Since K = kQ ■ Q'^, this implies that 

det C = (det Q)* = (c/2)a/A]A^ ivo 2 ql - “^vosq^qy + Voiq^), (48) 

where c is some unit amplitude complex number. This simple form implies that in the q —)■ 0 limit the zeroes of 
detC(q), which occur at [70] 


qy = 


qx 

Woi 


U03 ± VUo3 - ^^01^02 


(49) 
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depend only on Vg and not on the elastic moduli Ki and K 2 . The negative of the quantity under the square root is 
proportional to the determinant of the strain of the Guest mode: 

det£G = £xx£yy “ = Uq{voiVo2 - . (50) 

Thus, to order q^, there is a zero mode with real values of Qx and Qy if detec < 0 and complex values when deteg > 0. 
The complex values of q correspond to decaying or growing surface modes. The zeros at real values of q can either 
become complex when higher-order terms in the g-expansion are included, in which case they again correspond to 
surface modes albeit with an inverse decay length quadratic rather than linear in the surface wavenumber, or they can 
remain real. In the latter case, they can occur at specific values of q, in which case they are isolated. These points, 
which have a topologically protected winding number are called Weyl points, and they have received considerable 
attention recently in relation to Dirac semimetals and topological insulators [53, 71-78] and in photonic materials 
[79-81]. The distorted square lattices of figure 3(b) has isolated Weyl points [70]. It is also possible to have Weyl lines 
rather than isolated points. These occur in distorted versions of the 3fi pyrochlore lattice [82] and gyroid photonic 
crystals [80]. 


V. EXAMPLES OF PERIODIC LATTICES 

In this section, we will study three simple lattices: the square, kagome, and twisted kagome lattices depicted in 
figure 3. They provide specific examples of the phenomena discussed in the preceding section. Section V A and VB 
rely heavily of reference [45] and section V C on reference [83]. 


A. Square Lattice 

1. The nearest-neighbor lattice 

The periodic square lattice (figure 6 ) with only nearest neighbor bonds, along bi = a(l, 0) and 62 = a(0,1), and one 
site and two bonds per unit cell is the simplest example of a periodic Maxwell lattice. Many of its properties follow 
without calculation from the observations of the previous sections. It consists of straight filaments in orthogonal 
directions, each of which develops an SSS when placed under periodic boundary conditions. Thus an Nx x Ny lattice 
has Ns = Nx + Ny states of self stress and 

A^o = + Ny (51) 

zero modes, which, apart from the trivial translation modes, are infinitesimal floppy modes in which any row or 
column is rigidly displaced as shown in figure 6 (c). These rigid displacements form a basis for the null space of C, 
and their Fourier transforms do as well. If filaments parallel to the x-axis are rigidly displaced, then Qx = 0, and 
the Fourier transform of a set of displaced filaments parallel to the a;-axis are indexed by a wavevector <71 = {0,qy). 
Similarly rigidly displaced filaments parallel to the y-axis are indexed by q 2 = {qx,0)- There are Nx independent 
values of qx and Ny of qy. Each wavevector represents a zero mode, and since there are Nx independent values of 
qx and Ny independent values of qy, these account for all of the Nx + Ny zero modes, linear combinations of which 
account for rigid translations and rotations. Thus, the bulk phonon spectrum for the periodic square Maxwell lattice 
has two lines of zero modes running from the center F of the Brillouin Zone to the midpoint M and M' at the zone 
boundaries as shown in figure 7. 

An Nx X Ny lattice with free boundary conditions can be created by cutting and removing Nx + Ny bonds from 
an Nx X Ny lattice under PBCs. Thus, the cut lattice has only 2NxNy — {Nx + Ny) bonds [figure 6 ]. In the cutting 
process, however, all Nx -b Ny SSSs are lost, and the zero mode count of the free lattice is the same as that of the 
periodic lattice. To infinitesimal order in the displacement, the modes themselves are identical in the two cases. In 
the free lattice, however, the modes are nonlinear as shown in figure 6 (c). The bulk zero modes, which are seen under 
periodic boundary conditions, exhaust the Index theorem count. There are no additional zero modes at the surface 
of the free lattice. It is clear that any distortion of one of the surfaces parallel to the x- or y-axis of the free lattice 
will be transmitted across the entire sample by the rigid filaments, which support the SSS under periodic boundary 
conditions. So, zero modes are not localized near the surface; surface distortions have infinite penetration depth and 
thus do not constitute surface modes. 

There are exactly two SSSs at q = 0, and they clearly overlap affine bond elongations, which are equal to aSxx for 
bonds parallel to the x-axis and to aSyy for bonds parallel to the y axis. Thus the elastic energy density is simply 

fei = lk{elx+£ly). (52) 
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FIG. 6 : (a) An Nx = 5 by Ny = 5 square lattice with NN (full lines) and NNN (dashed lines) bonds showing bond directions 
for the calculation of self stress. The gray bonds at the right and upper boundaries are the N^ + Ny = 10 A and 2{Nx+Ny — 1) 
NNN bonds that must be cut to from a N^ x Ny lattice under periodic boundary conditions to produce the N^ x Ny lattice 
under free boundary conditions, (b) A unit cell showing its two Bravais lattice vectors ai and 02 and it two NN (1 and 2) 
and two NNN (3 and 4) bonds, (c) A finite mechanism of the finite square lattice with NN bonds only. This is also an 
infinitesimal mechanism in the periodic lattice. 


Qy 



M 

R 


r 

M' 





(a) 


FIG. 7: (a) Density plot of the lowest-frequency mode of the 
the Brillouin zone center at F to the midpoints M and M' of 
independent of qy and equal to zero for = 0. After reference 



NN square lattice showing lines of zero modes running from 
the zone edge, (b) 3d plot of the single mode ujx{q), which is 
[ 45 ] 


There are two independent lattice distortions, Sxx and Syy^ that cost energy and one Cxy that does not in agreement 
with the analysis of section IV B 1. 


2. The next-nearest-neighbor lattice 

Introducing NNN bonds, 63 = a(I, — 1 ) and 64 = a(l, I), increases the bond number without changing the number 
of sites. If NNN bonds are added one at a time, initially no additional SSSs are created, and each additional bond 
decreases the zero-mode count by one. As additional bonds are added, eventually additional states of self stress are 
created, for example in isolated cells with two NNN bonds in configurations like that of figure 4(d). Consider for 
simplicity the case with Nx = Ny. If a filament of equivalent contiguous NNN bonds (pointing either up or down) 
traverses the sample, a new SSS along that line is created. Thus, if there are no other NNN bonds, the change in 
bond number and number of SSSs relative to those in the state with no NNN bonds are, respectively, NNb = Nx 
and NNs = 1, leading to a decrease in the number of zero modes of AVg = —Nx -\- 1 for a total of Aq = Nx -\- 1. 
If each unit cell has an upward pointing NNN bond and not a downward pointing one, NNg = N^. In addition, 
there is one SSS for every q except q = 0 for which there are three for a total of Ns = N^ -b 2 leading to Nq = 2, 
i.e., leaving only the required two zero modes of uniform translation. Each added down pointing bond increases Nb 
and Ns by one leaving Nq fixed at two. Thus, if all NNN bonds are present, there are no zero modes beyond the 
two at q = 0, and the spectrum is gapped everywhere except at q = 0 as shown in figure 8(a). The gap at points M 
provides a characteristic frequency, 



oj* =2 


(53) 
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FIG. 8: (a) Comparison of the square-lattice phonon frequencies along symmetry directions in the BZ with fc' = 0 (dashed 
lines) and fc' > 0 (full lines) The blue curves depict ijJi{k',q), the lower and the red curves, uj 2 {k',q), the higher of the two 
modes. The single blue dashed line from i? to F represents the curves aJi(0,q'), uji{kl,q), and LjJ2{0,q), all of which are the 
same, (b) Frequencies for k' — 0 and fe' > 0 along the BZ edge from RtoM and back to R. When k' = 0, the frequency grows 
linearly with q away from M. When k' > 0, there is a gap defining the characteristic frequency w* = and a length scale 
when I* ~ Ijuj*. 


and a comparison of the fc' = 0 linear dispersion with the fc' yf 0 dispersion along the zone edge near M [figure 8 (b)] 
yields a characteristic length 


t 


a k 

2 V V' 


(54) 


which diverges as fc' —>■ 0. The approach of oj* and to zero as the Maxwell limit, fc' —)■ 0, is similar to the 

behavior of these quantities near jamming where oj* ~ (Az) and I* ~ (Az)”^ [37]. This result is not altogether 
surprising, given that weakening the spring constants oi NNN bonds to zero yields the Maxwell limit in which there 
is no characteristic frequency or inverse length. An effective medium analysis [46] of the case in which NNN bonds 
are randomly added with probability p = (z — Zc)( 4: leads to exactly the same results as jamming. 


3. Equilibrium and dynamical matrices 

We now investigate how these results follow from explicit calculations of the equilibrium and compatibility matrices. 
We set ai = o(l, 0) and 02 = a(0,1) and designate the NN bonds 1 and 2 to be to the right and above each lattice 
site, respectively, and NNN bonds 3 and 4 to be along ai — a 2 and ai -I- 02 , respectively, as shown in figure 6 (b). 
Following the rules outlined in section IIB, the force at site £ is f{Re) = ti{R() — ti{Ri — ai)-\-t 2 {Ri) — t 2 {Re — 0 , 2 ), 
and 


/I _ Q \ 

Q7VJv(q) = -(^ 0 (55) 

Thus if = 0 and qy ^ 0, the one-dimensional null space of Q(q) is spanned by the unit vector t^{qy) = (1,0); and 
if gy = 0 and qx ^ 0, it is spanned by the vector ty{qx) = (0,1). In other words, for each qy ^ 0, there is a SSS for 
each value of qy with independent tensions in bonds parallel to the a;-axis that have the same value for every bond in 
a given filament and similarly for qx yf 0. When both qx and qy = 0, there are two SSSs . Thus, there is one state SSS 
for each value of qx and one for each value of qy for a total of Ns = Nx + Ny SSSs . The null pace of C{q) = Qjvjv('7), 
which consists of the set of zero modes, is similar with one zero mode per qx and qy for a total of Nq = Ns zero modes 
as required by Index theorem. It consists of rigid displacements of individual rods as already discussed. 

The force at site £ arising from NNN bonds is f^^^{Re) = t^{Ri) — t-i{Ri — 63 ) -f t/i{Rt) — ti{Ri — ^ 4 ), and the 
NNN equilibrium matrix is 

1 / 1 — 1 _ p-i{q^+qv)u\ 

QNNN{q) — _|_ p-i{q^-qy)a ^ _ p-i{q:r+qy)aJ ' (56) 

When both NN and NNN bonds are present, there are four bonds per unit cell, and the full equilibrium matrix is 
a the 2x4 matrix 


Q{q) = {QNNiq) QNNNiq)) ■ 


(57) 
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FIG. 9: Densities of states piui) = p{uj)/p(Q) for the (a) square and (b) kagome lattices. The dashed black line (labeled 1) is 
the flat, one-dimensional pit^) at a; <C wo when k' = 0 , and the full red line (labeled 2 ) is p{uj) for k' > 0 showing linear-in-oj 
Debye behaviour at w <C w*, van Hove singularities near lo = uj* and constant behaviour at oj > cu*. Lines 3 and 4 in (a) 
are effective medium results for different probabilities of occupation of NNN bonds showing the washing-out of the van Hove 
singularity and the smooth transition from Debye to one-dimensional behavior. Adapted from references [45] and [46] 


At q = 0, all entries in Q are zero, and all q = 0 bond vectors are in its null space. Thus, the elastic energy is simply 
the expected affine result, 


4 

f — 2 V ^ ^ (58a) 

a—1 

= £yy) + 4 ^ [{Sxx + £yy + ‘^^xy) + {Sxx + £yy ~ “^^xy) ] (58b) 

= kifi^xx^yy + ‘ 2 K 44 sly, (58c) 


where Kn = k + k', and K 12 = K 44 = k'. As expected, shear moduli vanish when k' = 0. 
The NN and NNN dynamical matrices are 


D^^(q) 

D^^^(q) 


/sin^(q,ca/2) 0 \ 

m\ 0 sin^(qya/2)y ’ 

^k' A — cosqa;acosqyO sinsinq^a 
m I sinq^josinq^a 1 — cos qa;a cos qj,a 


(59a) 

(59b) 


The spectrum arising from the sum of these dynamical matrices is shown in figure 8 . When k' = 0, only 
contributes, and there are two independent one-dimensional modes with 


ujx,y = 2a;o| sin(q2,^j^a/2)|, (60) 

where wq = y^kjm. For every q^, u>x(q) goes to zero with as cQx, where c = wqu is the longitudinal (compressive) 
sound velocity, and reaches a maximum of 2a;o at the point M on the zone edge [q = (tt/o, 0)] and similarly for ujy{q) 
as shown in figure 8 . These one-dimensional spectra produce [65] (pages 205-209) a density of states 

1 / ^ \ 

p { u !) = - - ( 1 ~ —9 ) ) ( 61 ) 

27ra;oa \ dwg J 

which approaches a constant as w —)■ 0 and diverges as w —>■ 2uJo at the zone edge. 

When k' > 0, modes [with q = q(cos0,sin0)] exhibit a cos(40) angular modulation at low frequency and one¬ 
dimensional fc' = 0 behavior at larger q. When 0 < k' k, Dij{q) is well approximated as a diagonal matrix with 
mDxxiQ) = kq^ + Ak' s\t?{ qy/2) with associated eigenfrequency ojxiq) ~ \/Dxx{q)- These expressions immediately 
define the characteristic frequency of (53) at the point M = (0, tt) on the BZ edge. The first term in Dxx{q) represents 
the long-wavelength one-dimensional A^Wmodes that are present when k' = 0, whereas the second represents the 
effects of NNN coupling. When qx = 0, the only length scale in the problem is the unit lattice spacing, and no 
divergent length scale can be extracted from Dxx{0, Qy)- When the first term is large compared to the second, Dxxig) 
reduces to its form for the Maxwell fc' = 0 limit, and we can extract a length by comparing these two terms. The 
shortest length we can extract is that of (54), which comes from comparing kq^ to Dxx{q) at point M on the zone 
edge as depicted in figure 8 If q^ < tt, the Maxwell limit is reached when qx > q* ■ A similar analysis applies to Dyy{q) 
















17 



FIG. 10: (a) An = 5 by Ny = 5 kagome lattice showing NN (full lines) and NNN (dotted lines) bonds. The gray bonds 
along the right and upper edges are the 2{Nx. + Ny — 1) NN and 4:{Nx + Ny — 1) NNN bonds that must be cut from a Nx x Ny 
lattice under periodic boundary conditions to produce the the free Nx x Ny lattice, (b) Representation of a zero modes, (c) 
and (d) two different symmetric versions of the kagome unit cells showing labeling of sites and NN bonds. The vectors ai and 
02 , ai and — 03 , or any other similar pair can serve as basis vectors for the triangular Bravais lattice. 


when Qy > q* . If a square of length I is cut from the bulk, the wavenumbers of its excitations will be greater than 
ir/l, and for ql* > 1, all modes within the box will be effectively those of the lattice without NNN bonds. This 
construction is equivalent to the cutting argument of Wyart et al. [18, 38]. 

The characteristic length of (54) is identical to the length at which the frequency of the compressional mode 
0Jx{qx ~ 1/0 = Vk/l* y/Kii/l* becomes equal to uj*. A meaningful length from the transverse mode ijJx{0,qy) 
cannot be extracted in a similar fashion. The full phonon spectrum [figure 8(a)] exhibits acoustic phonons identical 
to those of a standard square lattice at g <C 1 and a saddle point at the point M. Thus, the low-frequency density of 
states shown in figure 9(a) is Debye-like: p{u}) ~ ujjy/kk' with a denominator that, because of the anisotropy of the 
square lattice, is proportional to the geometric mean of longitudinal and transverse sound velocities rather than to a 
single velocity. In addition p(w) exhibits a logarithmic van Hove singularity at uj* and approaches the one-dimensional 
limit {l/TT)/\/k at UJ* UJ 2y/k. The frequency uj* [(53)] is recovered by equating the low-frequency Debye form 
at UJ* to the high-frequency isostatic form of the density of states. 


B. Kagome lattice 

The NN kagome lattice consists of three grids of straight parallel filaments intersecting at lattice sites as shown in 
figure 10. This figure also shows two different unit cells reflecting the 3-fold symmetry of the lattice. For the moment, 
we focus on lattices with Nx = Ny cells on a side as shown in the figure. C derives the compatibility matrix for a 
generalized kagome lattice, of which the simple kagome lattice considered here is a special case. The equilibrium and 
dynamical matrices are straightforwardly calculated from it, as are the phonon spectrum and zero modes. As in the 
square lattice, each of the kagome-lattice filaments supports a SSS under periodic boundary conditions (care must 
be taken to join equivalent sites at the boundaries to create a single filament for bonds slanting to the left from the 
bottom), and the expectation is that a periodic lattice will have ‘tiNx states of self stress, and this is indeed the case. 
There is one state of self stress for each wavevector q = qGj /Gq along the symmetry equivalent lines from T to M in 
the BZ [figure 11(a)] parallel to the three reciprocal lattice vectors Gj, where Go = \Gj\ = 47r/(-\/3a). Since there are 
Nx values of q along each of these directions, there are a total of 3Nx SSSs, from which SSSs for individual filaments 
can be constructed. 

The 3Nj; SSSs require an equal number of zero modes, which, as in the square lattice, occur along lines in reciprocal 
space that have no component parallel to one of the three grids, i.e., along the lines TM in reciprocal space as shown 
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FIG. 11: (a) Shortest reciprocal lattice vectors, related by 3-fold rotations, of the kagome lattice that satisfy Gj ■ aj = 0. (b) 
Dispersion of lowest-frequency mode, showing “knife-edges” along the F — M in the BZ. (c) Density plot of lowest-frequency 
mode with k! jk = 0.02. Note the isotropic behavior near q — 0. Adapted from reference [45] 



(c) 



FIG. 12: (a) Phonon spectrum for the three lowest modes of the undistorted kagome lattice. Dashed lines depict frequencies at 
k' = 0 and full lines at k' > 0. The inset shows the Brillouin zone with symmetry points F, M, and K. (b) Phonon spectrum 
of the twisted kagome lattice with a > 0 and k' = 0. (c) Phonon dispersion along the zone edges from if to M in schematic 
form for both the undistorted lattice at fc' > 0 and for the twisted lattice with a ^ 0, showing the characteristic frequency 
UJ * = 2yF /m (uja ~ ^k/m\ sina|) and length Z* = (q*) ^ = (a/2)(fc'/fc)^/^ (Zq, ~ l/|sina|) for the untwisted (twisted) 
kagome lattice. Adapted from references [45] and [83[. 


in figure 11(b). The zero modes for a filament parallel to the a;-axis consist of displacements of all 1-sites and 3 sites 
by s(cos7r/6,sin7r/6) and s(cos7r/6, — sin7r/6), respectively, for infinitesimal s. This corresponds to rigid rotations 
of triangles about site 2 as shown in figure 10(b). An alternative description of the mode is that the entire filament 
is displaced a distance scos7r/6 to the right, and sites 1 and 3 are, respectively, displaced upward and downward a 
distance ssin7r/6 producing the zero-modes structure of figure 5(b). The spectrum of the lowest-frequency modes has 
a linear dispersion with uj = cq {c = y/ika/S) in the direction perpendicular to the T — M zero modes (figure 12(c)). 

As in the square lattice, adding NNN bonds gaps the spectrum leading to a characteristic frequency oj* ~ and 
associated length scale (* ^ Ij^/W calculated from the dispersion along the line M to K at the zone edge (figure 12). 
Other characteristic frequencies can be calculated from the lowest frequency optical modes or from the frequency at 
which the low-frequency acoustic phonon modes crosses over to a nearly flat dispersion. 

A total of ANx — 1 bonds must be cut to liberate a x Afr-unit-cell free lattice from its periodic parent. There are 
no states of self-stress in the free lattice, so there must be as many zero modes as bonds that are cut. This is more 
zero modes than the 3 A(e in the periodic system. The difference between the two numbers arises from the joining of 
lines slanting to the left under periodic conditions. The number of horizontal and right slanting filaments is the same 
in both cases. However, under PBCs, there are distinct left-slanting filaments; under free BCs, there are such 
lines terminating at the bottom and — 1 terminating at the right side of the lattice. 

Because of the three sets of filaments aligned along a„, at q = 0, there are now three rather than two SSSs, 
characterized by bond vectors ti = (1, 0, 0,1,0,0)/-\/2, t 2 = (0,1,0,0,1, 0)/-\/2, and ti = (0,0,1,0, 0, l)/-\/2, each 
of which has a nonzero overlap with the vector of affine distortions. This gives enough SSSs to fully stabilize the 
elastic energy of the NN kagome lattice with nonzero Lame coefficients [84]. Addition of NNN bonds increases these 
coefficients. 


8 


(fc -b 3fc') 


X = pi 


(62) 
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FIG. 13: (a) Twisted kagome lattice showing the displacements of sites in the two unit cells shown in figure 10 and rotation 
of triangles through ±q. Bond 1 (connecting sites 3 and 1) in this fignre is rotated through an angle —a. The lattice spacing 
for bonds of length a/2 is an = a cos a. (b) Lattices with different values of a superposed showing how changing a reduces the 
lattice area. From reference [83] 


The response is affine even though three sites per unit cell introduces the possibility of their undergoing nonafhne 
displacement to lower energy. However, the geometry of this lattice is special, and response is necessarily affine. 


C. Twisted kagome lattice 

The twisted kagome lattice is constructed from the finite zero modes by oppositely rotating triangles along all of 
the filaments of the untwisted lattice through an angle a as shown in figure 13. There are a continuum of lattices 
indexed by the angle a that bond 1 makes with the a;-axis. As we shall see, this lattice has properties that at first 
blush seem surprising but that in fact are simple consequences of the Index theorem . 


1. Bulk and elastic properties 

The straight filaments of the untwisted kagome lattice become “zigzagged” and lose their ability to sustain SSSs . 
The result is that there are only the two q = 0 states of self stress required by the Index theorem and the existence 
of two zero modes of translation, and there are no zero modes other than those at q = 0. Thus the simple rotation 
of triangles to create the twisted lattice from the the untwisted ones gaps all but the q = 0 bulk phonons, just as 
does adding NNN bonds to the untwisted lattice (figure 12). The untwisted spectrum is approached continuously as 
a —)■ 0 leading to a characteristic frequency (measured by the gap at the symmetry point M) and associated length 
scale. 


-- |sina|, L ^ -r-. -r- (63) 

I sin a I 

that, respectively, vanish and diverge as a —>■ 0. 

As figure 13 shows, twisting the lattice uniformly compresses it. If bond lengths are fixed at a/2, the Bravais 
lattice vectors are reduced in length from a to ol = a cos a, and the volume of each unit cell from (•\/3/2)o^ to 
(V3/2)a2 cos^ a. Thus angle changes modify the area of the lattice without changing any bond length of a NN 
lattice, implying that the bulk modulus B of these lattices vanishes for all a ^ 0. Observe that the twisted lattice has 
the peculiar property that it expands or contracts isotropically at no energy cost. If it expands in one direction, it 
will also do so in the opposite direction. Elastic materials with this property have a negative Poisson ratio [59]; they 
are auxetic. The twisted kagome lattice has the most negative allowed Poisson ratio of —1, which it retains for all 
strain, and it sometimes called maximally auxetic [83]. In addition, the unit cell contracts isotropically. Such lattices 
are termed equiauxetic [85]. There are not many naturally occurring materials that have this property, but artificial 
ones can be created [86]. 

As indicated above, there are two q = 0 states of self-stress that have the potential to create non-vanishing elastic 
moduli. The long-wavelength elasticity is necessarily isotropic, so if the bulk modulus is zero, the only option if for 
there to be an isotropic shear modulus or for none of the elastic moduli to be nonzero. The two states of self-stress 
overlap with affine strain, the shear modulus /r = ^/Sk/S is nonzero (and curiously independent of a and identical to 
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that of the untwisted kagome lattice [(62)], and the elastic energy density is 

/el “ 

where = Sij — ^Skk is the symmetric, traceless shear strain tensor. 

The strain is related to the metric tensor via Sij = {gij — dij)j2. The traceless part of the strain, iij = {l/2){gij — 
^Sij9kk), which is zero for gij = 6ij, is invariant, and thus remains equal to zero, under conformal transformations 
that take the metric tensor from its reference form Sij to h(x.)Sij for any continuous function ft.(x) of position x. The 
zero modes of the theory thus correspond simply to conformal transformations, which in two dimensions are best 
represented by the complex position and displacement variables z = x + iy and w{z) = Ux{z) + iUy{z). All conformal 
transformations are described by an analytic displacement field w(z). Since by Cauchy’s theorem, analytic functions in 
the interior of a domain are determined entirely by their values on the domain’s boundary (the “holographic” property 
[87]), the zero modes of a given sample are simply those analytic functions that satisfy its boundary conditions. For 
example, a disc with fixed edges (u = 0) has no zero modes because the only analytic function satisfying this FBC 
is the trivial one w(z) = 0; but a disc with free edges (stress and thus strain equal to zero) has one zero mode for 
each of the analytic functions w(z) = Unz'^ for integers n > 0. The boundary conditions limx^oo u{x,y) = 0 and 
u{x, y) = u{x+L, y) on a semi-infinite cylinder with axis along x are satisfied by the function w{z) = 
when = 2n7r/L, where n is an integer. This solution is identical to that for classical Rayleigh waves [59] on the 
same cylinder. Like the Rayleigh theory, the conformal theory puts no restriction on the value of n (or equivalently 
Qx). Both theories break down, however, at qx = Qc ^ niiii(^Qo~^) beyond which the full lattice theory, which yields 
a complex value of qy = q'y + iq'y, is needed. 


2. Surface modes 

As we have seen free two-dimensional lattices of N sites cut from a periodic Maxwell lattice necessarily have of 
order y/N zero modes because of order y/N bonds must be cut, and any sample-spanning states of self stress are 
lost under the cut. In the untwisted kagome lattice, these modes are identical to the bulk zero modes calculated 
under PBCs. In the twisted lattice, whose cut lattice must have the same number of zero modes as the untwisted 
lattice, there are no bulk zero modes (except at q = 0 ), and as a result the zero modes must be localized at surfaces. 
In the long-wavelength limit, these modes must reduce to the zero-frequency Rayleigh waves of an isotropic elastic 
continuum with vanishing bulk modulus with decay length Is = k~^ equal to the inverse surface wavenumber q. At 
shorter wavelength, Ig is determined by the length 1^ associated with the twisted phonon gap. 

Figure 14 depicts a finite rectangular sample with free horizontal surfaces parallel to ai and vertical surfaces 
parallel to 02 — 03 along with unit cells constructed so that all sites and bonds on a surface lie in periodically repeated 
continuous cells. It also shows which bonds (or bonds and sites) must be removed to liberate the finite lattice from 
the one under PBCs. The removal of two bonds or four bonds and one site per unit cell liberate the horizontal 
surfaces. In either case, the number of zero modes per cell is no = 2An — An^ = 2, where An = 0, An;, = —2 in the 
first case and An = — I, Anf, = —4 in the second. Similar arguments yield no = 4 for vertical surfaces. As might 
be expected, the modes are distributed equally between opposite surfaces, i.e., there is one zero mode per unit cell 
on horizontal surfaces and two modes per unit cell on vertical surfaces. Equivalently, there is one (two) mode per 
surface wavenumber — tt/Gs < q < tt/Gs, where Gg is the magnitude of the surface reciprocal lattice vector, 27 r/a for 
horizontal and 27r/(-\/3a) for vertical surfaces. 

The amplitude of the surface modes decay as exp[—its] with distance s away from the surface all of the way to the 
opposite free surface, where k is in general complex indicating oscillations along with decay, k = Rek = is the 
inverse decay length. In the case of horizontal surfaces, the decay length is the same on opposite surfaces. Figure 
14(b) plots the single «:(q, a) for different values of a. In the case of vertical surfaces, the two decay lengths for the 
left surface differ from those of the right one. Figure 14(c) plots these for a = tt/S. In all cases, one «;((?, a) reduces 
to K{q, a) = q in the long-wavelength limit as required by the continuum theory. 

Surface zero modes are by definition in the null space of the compatibility matrix C. Systems with parallel free 
surfaces with PBCs along their direction of alignment can be viewed as a series of layers n = 1, and C can be 

decomposed as 
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FIG. 14: (color online) (a) A free twisted kagome lattice with free horizontal and vertical boundaries. Sites and bonds of unit 
cells that match the surface are shown in black and red, respectively, and the geometric form of these unit cells are indicated 
by dashed quadrilaterals. Cutting the two (four) dashed bonds per cell on the horizontal (vertical) boundary produces a lattice 
with bottom (left) and top (right) boundary cells 1 (4) and 2 (5). Alternatively removing these bonds along with one circled 
site and the two bonds crossed by grey lines changes the top (right) boundary cell from 2 (5) to 3 (6). The number of zero 
modes per site is 1 (2) per site on both the bottom (right) and top (left) boundaries. Cell 7 has one additional bond cut from 
it. (b) The reduced inverse decay length Kr = k/(Gt./2) of the horizontal boundaries as a function of qr = q/{Grl2), where 
Gr = 27r/a for a = 7r/20, rr/lO, 37r/20, rr/S, 7r/4 in order from bottom to top. All curves follow Kr — qr near qr = 0. The curve 
at a = 7r/4 diverges at qr = 1.0. (c) Kr as a function of qr with Gr = 27r/\/3 for a = tt/S. The positive curves are for the left 
boundary and the negative ones are minus Kr for the right boundary. The two positive curves merge at qr « 0.4. There are 
still two distinct decays beyond this point with the same real part and different imaginary parts. The straight grey lines are 
the elastic limit Kr = qr- 


where all C^ ’s depend on the surface wavenumber q and a. Cn is the 6x6 matrix connecting bonds and sites within 
in a single unit cell and C 12 is the 6x6 matrix connecting bonds in one unit cell to sites in unit cells one layer deeper 
in the sample. The same unit cells are used throughout the sample. The opposite surface may terminate with only 
a partial version of these cells, and the exact forms of and Cnn depend on how the surface is terminated. 

Consider, for example, horizontal surfaces, the bottom of which is characterized by the unit cell labeled 1 in figure 
14(a) and the top of which is characterized by unit cell 2. For modes localized at the bottom surface, unit cell 1 is 
used through out the sample. The termination cell 2 at the opposite surface has all three sites but only four bonds of 
unit cell 1. The missing bonds (5 and 6) are not affected by cell N — 1. Therefore, Cjv,Ar-i = C 12 is a 6 x 6 matrix, 
and CATTvisadxb matrix. If the top surface is terminated by unit cell 3, which has only 2 sites and 2 bonds of cell 
1, Cjv-i,Ar is a 6 X 4 matrix and Cn,n is a 2 x 4 matrix. Thus displacements U = {ui, ...,m^) in the null space of C 
satisfy 


CiiUn + Ci2M„+i — 0, (66) 

for n = 1, ...Af — 2. These equations are solved by = Am„ and 

det(Cii + AC 12 ) = 0 (67) 

subject to the boundary conditions that CuMat-i + Cn-i^nUn = 0 and = 0. 

The inverse penetration depth is determined by A via exp[—«;aj_] = A, where aj_ = a-\/3/2 is the distance between 
unit cells in the direction perpendicular to the surface. In case of termination with either unit cell 2 or unit cell 3, 
un must equal Amjv-i to solve the first boundary condition. Though required by the known existence of the zero 
mode, it is a remarkable fact that the projection of Ui, which is in the null space of Cn + AC 12 , onto the space of 
displacements of the last layer is also in the null space of Catat, and as advertised earlier, the surface mode decays 
exponentially from one free surface to the next. A similar analysis applies to the vertical surface for which a_L = a/2. 
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FIG. 15: (a) The zigzagged square lattice with two sites per unit cell and rectangular unit cells and Brillouin zone, (b) and (c) 
density plots of the two lowest-freqnency bands, showing lines of zero modes along qy = 0 and = 0 in (a) and along = 0 
in (b). The total nnmber of zero modes in the two bands is the same as in the undistorted sqnare lattice. 


Of course, any linear combination of the exponentially decaying modes on the two sides of a strip is a zero mode if 
both are individually. In the usual situation in which the Rayleigh waves have a nonzero frequency, the eigenstates 
of a finite strip are symmetric and anti-symmetric combinations of states on the two surfaces that interact across 
the strip yielding a peristaltic mode with a; ~ g and a bending mode with w ~ It should be noted that k for 
modes that are not localized within the surface unit cell can also be calculated from detC'(< 7 ), which is a function of 
exp[iq • ai], i = 1,..., 3, by setting q = q±'hin + q\\ where < 7 || is the component of q parallel to the surface and fiin is 
the unit inward normal to the surface, setting q± = ik, and solving for k in detC(jK, q||) = 0. This approach does 
not directly provide eigenvectors satisfying boundary conditions. 

A surface consisting of unit cells 7 in figure 14(a) differs from surfaces composed of other unit cells shown in that 
figure in that it is missing a bond on the surface: it is obtained from cell 1 by cutting the dashed downward pointing 
bond, and as a result, this surface has two, not one zero mode per cell. The calculation just outlined indeed produces 
two zero modes per q, one of which is localized completely on the first row because Cn has a non-empty null space. 
Similarly, k diverges as shown in figure 14(b) at q = tt/Gs for a = 7r/4 because Cii(7r/Gs, 7r/4) has a non-empty null 
space. 


D. Other lattices 

So far, we have focussed on the three simplest examples of two-dimensional periodic Maxwell lattices and the 
free lattices cut from them. There are many others that can be constructed from these without changing the local 
coordination numbers or the lengths of any bonds and whose properties can easily be understood in the context of the 
Index theorem. One of the simplest of such lattices is the “zigzagged” square lattice with two sites per unit cell [88], 
shown in figure 15, in which every other row is displaced to the right while allowing the requisite compression along 
the vertical direction. This lattice retains the straight horizontal filaments of the original square lattices but loses 
those in the vertical direction. It does not, however, lose any vertical SSSs because the SSSs from individual straight 
filaments are converted to ones like those of figure 5(e) on pairs of zigzagged filaments of which there are total of 
under PBCs . These SSSs must be accompanied by zero modes in the spectrum, which show up as horizontal and 
vertical lines of zeros in the lowest-frequency mode and a vertical line of zeros in the second-lowest frequency mode 
as shown in figures 15(b) and (c). The vertical SSSs have no overlap with affine strain, and the lattice does not resist 
vertical compression. The result is that the elastic energy density is simply /ei = ku^^/2. 

The kagome lattice offers more interesting variations [83, 88]. Figure 16(a) shows one such variation with intriguing 
properties. It consists of alternating rows of oppositely tilted distorted hexagons. It has rectangular symmetry with 
6 sites per unit cell. As figure 16(b) shows, it has an unusual spectrum: its modes are fully gapped (except at q = 0) 
near the origin but exhibit curved lines of zero modes at large q. It has zero-frequency surface modes for surface 
wavenumbers in the gapped region but not in the ungapped region. Curiously, the elastic energy density is identical 
to that of the twisted kagome lattice even though it has a lower symmetry. 

The square and kagome lattice are the two-dimensional Maxwell lattices with the smallest unit cells. As discussed 
above, the kagome lattice can be distorted in various ways to produce larger-unit-cell Maxwell lattices, but there are 
many other lattices including random ones [89] that can be created. An intriguing set of periodic Maxwell lattices [69] 
are those arising from rational approximates to Penrose rhombohedral tilings [90] that approach a quasi-crystalline 
lattice with 5-fold symmetry [91-93]. A unit cell for the second periodic approximate with 80 sites and 160 bonds per 
cell is shown in figure 17. These periodic lattices, which can be constructed via projection from a five-dimensional 
cubic lattice [94], all have an average coordination of exactly four even though the coordination of local sites varies 
from three to as high as ten. The size of the unit cell increases rapidly from = 30 in the first approximate to 
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FIG. 16: (a) An example of one of the many lattices that can be constructed from the kagome lattice without changing any 
bond lengths overlayed on an undistorted kagome lattice. This one is a stacked lattice with alternating rows of oppositely tilted 
hexagons, (b) Density map of the lowest-frequency band showing the gapped spectrum near the origin and two curved lines of 
zero modes passing between two corners. Adapted from reference [83] 



FIG. 17: (color online) Unit cell of the second periodic approximant of the Penrose tiling showing SSSs (circled 1, 2, and 3). 
In all states, stress is localized on vertical ladders with different signs of stress on opposite sides as in figure 5(c). States 2 and 
3 share the bond marked bs, and are not orthogonal. They can be orthogonalized to produce states mostly, but not completely 
localized on the two ladders. From reference [69] 


N = 25840 in the eighth approximate. Each of the approximates is a legitimate periodic lattices with a full set of 
phonon branches with dispersions depending on lattice wavenumber q. They can also be interpreted as a single-cell 
system under toroidal PBCs that approach the infinite-cell quasicrystalline limit. In this interpretation, which we 
pursue here, only the q — 0 part of the spectrum is of physical interest as is the case for periodic approximates to 
randomly packed spheres at jamming. These lattices have a number of interesting properties: (1) their undistorted 
versions have of order y/N SSSs and zero modes, but none of the SSSs overlap with affine strain, and all elastic 
moduli are zero, much as in rigidity percolation. (2) Randomizing site positions removes all but the two required 
SSSs and zero modes. The two SSSs overlap with affine strain to produce two independent positive eigenvalues of 
the elasticity matrix, one corresponding in the large N limit to the bulk modulus with associated eigenvector of pure 
compression and one to a shear modulus with associated shear eigenvector. The bulk modulus increases with N 
reaching a saturation value at N 10'^. The shear modulus on the other hand approaches zero as 1/7V. The latter 
behavior is required by 5-fold-symmetric quasicrystalline limit whose elasticity must be isotropic with both shear 
moduli equal. Because there are only two SSSs, one shear modulus must be zero if the bulk modulus is nonzero, and 
the other must approach zero with N. This is essentially the same behavior observed in randomly packed spheres at 
the jamming threshold with z = 2d [34, 35]. 


VI. TOPOLOGICAL PHONONS 

Twisting the kagome lattice gaps all zero modes of the untwisted lattice except those at q = 0. This gapping is 
reminiscent of that of the electron spectrum in systems like polyacetylene [52, 95], quantum Hall materials [50, 51] 
and topological insulators [53-58], which is associated with the appearance of topologically protected zero modes at 
free boundaries and at boundaries separating two topological classes. An interesting and natural question is then 
whether or not the phonon spectra of Maxwell lattices can be gapped in different ways to produce distinct topological 
classes with protected boundary modes. The answer is yes they can be, and this section, which is mostly based on 
reference [60], will explore both how they can be constructed and the nature of their interface states. 
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FIG. 18: (a) and (b) depict the SSH model of polyacetylene, with A and B sublattices indicated by blue circles and red sqnares, 
respectively, (a) describes the gapless state with all bonds identical, while (b) describes the gapped AB dimerized state, with 
donble (single) bonds on the AB (BA) bonds. The BA dimerized state with single and double bonds interchanged is not shown, 
(c) and (d) show the ID isostatic lattice model in which masses, represented by the larger blue dots, are connected by springs 
in red and are constrained to rotate about fixed pivot points, represented by small black dots, (c) is the gapless high-symmetry 
state with 0 = 0, and (d) is the gapped lower-symmetry phase with 6 > 0. (c) and (d) are equivalent to (a) and (b) if we 
identify the masses (springs) with the A (B) sublattice sites, (e) shows a domain wall in polyacetylene connecting the AB and 
BA dimerized states. There is a topologically protected zero-energy state associated with the A sublattice at the defect, (f) 
shows the equivalent state for the mechanical model with a topologically protected zero-frequency mode at the domain wall 
connecting a. 9 = +9c lattice with a 9 = —9c lattice, (g) shows a domain wall connecting the BA and AB dimerized states, 
which has a zero energy state associated with the B sublattice, (h) shows the equivalent isostatic state with a state of self-stress 
(SS) at the domain wall and zero modes at the end so that Index count No — Ns = 1 is satisfied. Adapted from reference [60] 


A. A one-dimensional Model 

We begin with a one-dimensional model whose phonon spectrum (including both positive and negative frequencies) 
is identical to that of the one-dimensional Su-Schrieffer-Heeger model for polyacetylene [52, 95] schematically depicted 
in figures 18(a) and (b). Our model, depicted in figures figure 18(c) and (d), consists of rigid bars of length r that can 
rotate freely about fixed positions on a one-dimensional periodic lattice. The ends of neighboring bars are connected 
by harmonic springs whose lengths are adjusted so that the equilibrium configuration is one in which alternate rods 
make an angle 0 with the upward or downward normals. Bars tilt towards the right if 0 > 0 and to the left if 0 < 0. 
Each rod s has one degree of freedom 9^ = 9 — 56s, and each spring provides one constraint. Under periodic boundary 
conditions, the number N of rods equals the number Nb of springs. In the linearized limit, the compatibility matrix 
with components C^s connects the stretch in spring /3 with rotations of rod s: 5lp = CjSsSOg, and the equilibrium 
matrix Qs/s relates torques on rod s to tensions tp in spring /3: Ts = —Qsjitp- With appropriate sign convention for 
the torque, Cps = Q'^s- ^ state of self stress, springs are under tension, but there are no torques on the rods. The 

Index theorem thus applies directly to this system. Under periodic boundary conditions Nb = N and Nq = Ns- 
Cutting one bond creates a free lattice with Nq = Ns + 1. Thus, if there are no SSSs, the free system must have one 
zero mode, which can either be a mode in the bulk spectrum or a surface mode on one of the boundaries. But what 
determines whether it is localized on the right or left boundary? As we shall now, it is the topological class. 

We proceed now to a more detailed analysis of the our model. The components of the compatibility matrix at rest 
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angle 0 are 

Csp{6) = Ci{9)6a,p - C2(0)(5/3,s+1, 

where 

(a ± 2 r sin 9)r cos 6 
^ ^ \/a'^ + 4r2 cos^ 6 

Thus |ci| > |c 2 | for all 0 < 0 < tt, and |ci| < |c 2 | for all —tt < 0 < 0. The energy of the system (contained entirely in 
the stretching of the springs) is then 

E = = lkJ^(ciS0a - C2S0a+l f. (70) 

/3 

The Fourier transform of Cs /3 is 

C(q) = Cl - e“i"c 2 , (71) 

and bulk phonon modes have frequency 

w(g) = ±|C'(g)| = (ci - C 2)2 + 4 ciC 2 sin^(ga/2) (72) 

(for unit mass). When 0 = 0 (vertical rods), Ci = C 2 , the energy becomes invariant with respect of 60^ —> 60^ + 6 for 
every s, and there is necessarily a bulk zero mode at g = 0 - this in spite of the fact that the bases of the rods are 
anchored, breaking translational invariance. For other values of 0, the phonon spectrum is fully gapped. 

In a finite system, the compatibility matrix can be expressed in the form of (65) with Cn = ci and C 12 = —C 2 . 
The decay length of the surface state is determined by Ci — Ac 2 = 0 or A = Ci/c 2 = (e'^®) for states localized 

on the left (right) boundary. Thus, the one zero mode is on the left if |ci| < |c 2 | and on the right for |ci| > |c 2 |. 
The left zero mode is particularly easy to see at the critical angle 9c = — sin“^[a/(2r)] at which ci = 0 as shown 
in figure 18(s). The equilibrium matrix Q = for a finite system has no zero modes, and there are no SSSs as 
required by the Index count. This follows because Cuti = 0 for any t in the null space of Q. Then the equation 
—Ci 2 tn + Ciit„+i = 0 sets all for n > 1 equal to zero. Under periodic boundary conditions, there must be one 
localized state of self stress for each localize zero mode. 

But what does this have to do with topology? The compatibility matrix C{q) = \C{q)\'^’^ (or more generally its 
determinant) maps points in the Brillouin zone (—tt/o < <7 < tt/q) to a path in the complex plane. Since it depends 
on the path will be a closed and return to its starting point as q advances between equivalent points in the zone 
[q ^ q + (27r/a)]. These curves are characterized by an integer winding number: 

1 /• 2 ’r 1 cl-Kja 1 

n= — / d(j) =-— / dg—Imlndet < 7 ( 5 ), (73) 

27r dp 27ri Jq dq 

which for the C{q) of (71) is either +1 or 0. Clearly n = 0 if the path in the complex plane does not enclose the origin, 
and n = 1 if it does. The first case occurs whenever |ci| > |c 2 | and the second whenever |ci| < |c 2 | as shown in figure 
19. When |ci| = |c 2 |, the boundary of the curve passes through the origin. The winding number is thus a topological 
invariant equal to 1 for all —tt < 0 < 0 and to 0 for all 0 < d < tt . The only way it can change values is by passing 
through the critical angles 0 = 0 or 0 = tt (bars lie along the horizontal axis). When n = 0, the zero mode is on the 
right and if n = 1, the zero modes is on the left. The connection between the topological invariant and the existence 
of zero modes is easy to see in this case. A zero surface mode exists if there is a solution to c(A) = ci — C 2 A = 0 with 
|A| < 1. The compatibility matrix is simply C(q) = c(A = e^^*^), and along the closed path it describes in the complex 
plane, |A| = 1. In the interior of the path |A| < 1. Thus, if the path encloses the origin, which is the point at which 
c(A) = 0, the solution for A will have a magnitude less than one. A complementary perspective, based on the Cauchy 
argument principle, is discussed in D, 

Here we considered the winding number associated with C(q). We could equally well have considered that associated 
with the equilibrium matrix Q(q) = C*(q). Because winds clockwise rather than counter clockwise with 

increasing a, its winding number is either —1 or 0, with —1 corresponding to |ci| < |c 2 |. Thus the surface state 
is at the left when the winding number of Q is —I. We will use a generalization of the Q-winding number in our 
characterization of topological classes of kagome lattices in what follows. 

As is well known for the SSH model [52, 96], an interface between the two dimerizations binds a zero mode, as 
indicated. Similarly, the interface separating the two signs of 0 does as well. This is most easily seen for 0 = ±9c 


( 68 ) 

(69) 
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FIG. 19: Contour plots in the complex plane of C{q) for complete circuit of q from q — 0 to q = 2-71 /a. (a) C2 = 1.0 and 
Cl = 0.5 < C 2 ; the circuit contains the origin and the topological charge is n = 1. (b) ci = C 2 = 1.0; this is the transition case 
in which the contour just touches the origin, and n = 0. (c) Ci = 1.5 > C2 = 1.0; the circuit does not contain the origin and 
n = 0. Note the circuits begin with g = 0 on the far left and circulate counterclockwise. 


where the springs are colinear with the bars, so that ci or ci =■ 0. Figure 18(f) shows a domain wall between -\-9c 
and —0c, in which the center two sites share a localized zero mode. Figure 18(h) shows an interface between —0c and 
+0C with a state of self-stress localized to the middle three bonds, in addition to floppy modes localized at either end. 
As long as there is a bulk gap, the zero modes cannot disappear when 0 deviates from ±0c. The zero modes remain 
exponentially localized, with a localization length, I = a/ln|ci/c 2 |, that diverges when 0 —> 0. 


B. Topological lattices 

We have just seen the intimate connection between topological properties of the compatibility or equilibrium 
matrices and zero-modes at boundaries in a one-dimensional system. Here we discuss how these ideas can be extended 
to higher-dimensional Maxwell lattices. We will for the most part only quote results, and not provide details of how 
they were obtained. The latter can be found in reference [60]. There are two significant differences between these 
lattices and those in one dimension. First, because Q(q) (or C(q)) depends on the vector wavenumber q with two 
rather than one independent component, topological characterization will require a vector rather than a scalar winding 
number. Second there are boundary surface modes imposed by the Index theorem that are required whether or not 
lattices have any topologically properties. We will thus need to divide the surface mode count into those parts imposed 
by the Index theorem and those parts implied by topological considerations. The result is that topology cannot change 
the global count of the Index theorem, but it can move zero modes from one boundary to another and give rise to 
topologically protected zero modes at interfaces between two different topological classes. 


1. Topological and total mode count 

To unify the treatment of zero modes arising from mismatch of sites and bonds and those that arise in locations 
where there is no local mismatch, reference [60] generalized the Index theorem so that it determines a zero mode 
count 10 ^ = Nq — N§ in a subsystem S' of a larger system. This is well defined provided the boundary of S is deep 
in a gapped phase where zero modes are absent. This count has the two separate contributions just discussed: 

+ Vp, (74) 

where is a local count of sites and bonds in S and is a topological count, which depends on the topological 

structure of the gapped phases in the boundary of S. The topological contribution has a form similar to that of the 

one-dimensional system. For the periodic lattices we are considering, it is best expressed as a count per unit cell on 
an edge indexed by a reciprocal lattice vector G (i.e., G is normal to the surface with a magnitude of 27r/o_L, where 
a_L is the distance between lines of cells identical to the line of surface cells), 

i>^ = u^/N,,ii = G-RT/{2T:), (75) 

where A^ceii is the number of surface unit cells and Rt, a generalization of the one-dimensional winding number, is a 
lattice vector 


R’T — ^ ^ 


(76) 
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where are the primitive translation vectors and 

where (j){q) is the phase of detQ(< 7 ) {Q{q) = Here Ci is a cycle of the BZ connecting q and q + Bi, 

where B^ is a primitive reciprocal vector satisfying • Bj = 2TT6ij {Bi = —G 2 and B 2 = Gi in figure 11a). The 
rii are winding numbers of the phase of det Q(<7) around the cycles of the BZ, where Q(q) is the equilibrium matrix 
in a Bloch basis. This winding number is independent of path, and thus independent of q so long as the spectrum 
is gapped everywhere except the origin. The zero-mode at the origin of q is not topologically protected (i.e, it can 
be gapped by a weak potential breaking translational symmetry), so it does not cause any problems. It is possible, 
however, for there to be topologically protected gapless points. These would be point zeros around which the phase 
of det Q(k) advances by 27r. These lead to topologically protected bulk modes that form the analog of a “Dirac 
semimetal” in electronic systems like graphene [53, 71-78]. These singularities could be of interest, but they do not 
occur in lattices derived from the kagome lattice we study below. They do, however, occur in modified square lattices 
[70] of the type shown in figure 3 considered in reference [49] and in the pyrocholore lattice [82]. 

In general, the winding number is not gauge invariant and depends on how the sites and bonds are assigned to unit 
cells. It is, however, possible to adopt a “standard unit cell” with basis vectors T’^(,g) for the Ug sites (nt, = dus bond 
centers) per cell for which the “dipole” moment of the site and bond charges, i?stan = is equal to 

zero. The two unit cells of figure 10 satisfy this zero dipole condition even after being distorted to those of the twisted 
lattice. Q(q) is defined using Bloch basis states |q, a = n, (3) cx ^^expiq • (Ri -|- ra)|Ri -I- ra), where Ri is a Bravais 
lattice vector. In this gauge, Rt is uniquely defined, and the zero-mode count is given by equations (74) to (76). 

Because Rstan is zero, we could equally well use a basis in which all and rp are simply zero. This is in fact 
the basis we use for all of the calculations presented here. It should be noted that it is not always possible to find a 
symmetric “standard” unit cell with a vanishing dipole moment defined in terms of charges at sites and bond centers. 
Indeed, there is no such cell for the SSH model or for 3D pyrochlore lattices. Not having such a standard cell is not 
a problem, however. The number and the nature of surface modes do not depend on the choice of a “standard” or 
reference unit cell and are unambiguous. The easiest choice is usually to set the positions of all of the sites and bonds 
in the unit cell to be zero. As discussed in D, det C for different unit cell choices, such as those used to calculate 
the zero surface modes in section VC2, will have different phase factors, which account for the differences in the 
topological integral of (77) for different choices of unit cells.. 

The local count, depends on the details of the termination at the surface and can be determined by evaluating 
the macroscopic “surface charge” that arises when charges +d (—1) are placed on the sites (bonds) in a manner 
analogous to the “pebble game” [11]. This can be found by defining a bulk unit cell with basis vectors fa that 
accommodates the surface with no leftover sites or bonds (see figures 14 and 20) as discussed in section V C 2. This 
unit cell depends on the surface termination and, in general, will be different from the “standard” unit cell (figure 10) 
used for the calculation of The local count is then the surface polarization charge given by the dipole moment 
Rl per unit cell: 

H = lyf/Naeii = G ■ RL/2n, (78) 

where 

Rl = d 

sites jj, bonds 0 

The total zero-mode count on the surface then follows from equations (74), (75), and (78). The polarization of the 
standard cell is zero so that Rl = Rl — Rstan can be calculated from the displacements of the surface cell 

sites and bonds relative to those of the standard cell. 


2. Constructing topological lattices 

The kagome lattice has three sites per unit cell, and displacing these sites while maintaining the size and shape of 
the unit cell creates different lattices, each of which can be smoothly connected with the other across domain walls. 
The twisted kagome lattice, with bond length increased by 1/ cos a to connect smoothly with the untwisted lattice, 
is an example of this operation. The most general such lattice is described by four parameters (one site in the the 
unit cell can always be fixed). The most useful parametrization is one in which the “straightness” of the three sets 
of filaments are controlled individually. Such a set is depicted in figure 20. Displacing site 1 of the unit cell by 
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FIG. 20: (a) and (b) depict the geometry used to derive equations (80) and (81). In (a), site 1 of the symmetric unit cell of 
figure 10 (c) is displaced downward by ^/hxzpi, perpendicular to the lattice vector ai, and site 2 is displaced along as by 0 : 103 . 
The result is that the filaments along 02 and 03 remain straight whereas those along oi are zigzagged, (b) depicts similar 
displacements that zigzag only filaments along 02 or 03 . (c) shows the ternary phase diagram for fixed oii + 2:2 + 0:3 > 0. The 
space enclosed by the central triangle corresponds to the to the state with Rt = 0. The direction of Rt is the six sectors 
surrounding the central triangle are indicated by the blue arrows. The red dots from right to left correspond to the three 
lattices shown in the opposite order in figure 21 


—where pi is the vector of length a perpendicular to lattice vector ai, “zigzags” the filaments parallel to 
ai- Though this operation leaves filaments parallel of 03 straight, it zigzags the filaments parallel to 02 ■ The latter 
can be straightened by the simple operation of displacing site 2 along 03 by xia^, as shown in figure 20(a). This 
process is repeated to zigzag filaments parallel of 02 and 03 as shown in figures 20 (b) and (c). Finally the 1 — 2 — 3 
triangle can be isotropically expanded by displacing the three sites the same amount along directions pi, P 2 and p^. 
The basis vectors for the unit cell are then 


rt, = r°- VSx^Pf, + cc^-ia^+i + {z/V3)Pt,-i, (80) 

where r° (e.g., Vi = ai/ 2 , r 2 = —a 3 l 2 , = 0 ) are the basis vectors of the untwisted unit cell and it is understood 

that all subscripts are evaluated mod 3. The bond vectors are then 

bf3 = dp+i -dp = \ap - (-l)^"h(/3-i)/3)[('3.^_^^ _ y/Zx^pp], (81) 

where Int(a;) is the integer part of x and it is understood that a/ 3+3 = 0/3 and P/ 3+3 = pp- The expressions for 
6/3 for /3 = 4,5,6 are obtained from those for /3 = 1,2,3 via the relation 6/3+3 = ap — bis. The untwisted kagome 
lattice corresponds to X = {xi,X2,X3;z) = (0,0,0;0) and the twisted kagome with twist angle a = tan“^(2-\/3a:) 
to X = (x,x,x; 0 ). For a lattice with straight filaments along Oi only, X = ( 0 , X2, X3; z) and similarly for straight 
filaments along 02 and 03 . These lattices have SSSs and associated zero modes along a single direction in the Brillouin 
Zone. Moving off a/i =0 gaps the spectrum as shown in figure 21. 

The topological polarization in terms of X is 


Rt 


1 X ^ 

-^ap sign(a:p). 

p^l 


(82) 


This leads to the ternary phase diagram shown in figure 20(c) as a function of xi, X 2 , and X 3 for fixed xi + X 2 + X 3 
and z = 0. It has eight octants corresponding to the eight possible sign combinations of (xi,X 2 ,X 3 ). The (+,+,+) 
and (—, —, —) octants correspond to the class of the twisted kagome lattice. The remaining 6 octants are states that 
are topologically equivalent, but related to each other by Cg rotations. Figure 21 shows representative structures for 
the Rt = 0 phase (figure 21 (a)), the Rt yf 0 phase (figure 21 (c)), and the transition between them (figure 21 (b)). 
The insets show density plots of the lowest frequency mode, which highlight the gapless point due to the acoustic 
mode in figure 21(a) and the gapless line due to states of self stress in figure 21(b). In figure 21(c), the gap vanishes 
only at the origin, but there is a low-frequency cross that arises because acoustic modes vary quadratically rather 
than linearly with q along its axes. This behavior will be discussed in the next section. 

We are now in a position to analyze zero interface modes for different X. Surfaces in the twisted-kagome octants 
are in the same class of those discussed in section V C 2 and in figure 14. The bonds and sites to be displaced to 










29 



FIG. 21: Representations of lattices and the density maps of their associated lowest frequency modes for (a) a twisted kagome 
lattice with X = (0.1,0.1,0.1; 0) and Rt = 0 (right-most red dot in figure 20(c)), (b) a critical lattice with X = (0, 0.1, 0.1; 0) 
(central red dot), and (c) a topological lattice with X = (—0.1, 0.1,0.1; 0) and Rt = —ai (left-most red dot). Note the 
isotropically gapped spectrum (except for the origin) of (a), the line of zero modes in (b), and the soft-mode cross in (c). 
Adapted from reference [60] 








FIG. 22 : Real part of the reduces inverse penetration depth Kr = 2kIGs for various orientations of surfaces and Rt as a 
fnnction of reduced surface wavenumber qr = 2 qlGs, where Gs is the magnitude of the surface reciprocal lattice vector. In (a) 
the surface lattice vector G = G2 — G3, Rt ~ —ai, G ■ RT/{2-n) = 2 , and there are four zero surface modes. The opposite 
surface with —G has no zero modes. The full magenta curve is doubly degenerate with opposite imaginary parts, and the 
two dashed curves are not degenerate, (b) G = — Gi, Rt = <13 and G • Rt/{2tt) = 1 . There are two surface zero modes, 
(c) G = G2 — G3, Rt = —0,2, G ■ Rt/{2tt) — 1 , and there are three surface zero modes. The full magenta curves in (a) to 
(c) correspond to acoustic surface states, and in all cases Kr approaches 0 as q] as qr —>■ 0. (d) shows how to construct the 
dipole moments Rl for the surface unit cells in (a) to (c). On the left, black bonds 4 and 6 of the symmetric unit cell are, 
respectively, displaced by 02 and —03 while the circled site 1 is displaced by 02 to produce the surface unit cell of (a) and (c) 
with Rl = 202 — {02 — 03) = —oi and Ul = 2. On the right, the black bond 6 is displaced through —03 to produce the surface 
unit cell of (b) with Rl = 03 ad vl = 1 . 


TABLE I: Reciprocal lattice vectors G indexing surfaces, dipole moment Rl, and local count ul for the seven surfaces cells of 
figure 14(a). The reciprocal lattice vectors Gi, G 2 , and G 3 are depicted in hgure 11 and the Bravais lattice vectors ai, 02 , 
and 03 in figure 10 . 


cell 

1 

2 

3 

4 

5 

6 

7 

G 

-Gi 

Gi 

Gi 

G 2 -G 3 

G 3 -G 2 

G 3 -G 2 

-Gi 

Rl 

03 

02 

-03 

-Oi 

Ol 

Ol 

03 — 02 

i>L 

1 

1 
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FIG. 23 : Zero modes at domain walls, (a) shows a lattice with periodic boundary conditions and two domain walls, the left one 
between (.1,.1,.1;0) and (.1, .1, —.1; 0) with zero modes and the right one between (.1, .1, —.1; 0) and (.1, .1,+.1; 0) with states 
of self stress. The zero mode eigenvectors at Qj, = tt are indicated for the floppy mode (arrows) and the state of self-stress (red 
(-f) and green (-) thickened bonds), (b) shows the vibrational mode dispersion as a function of q^- From reference [ 60 ] 


convert the standard unit cell to a surface one are depicted for two cells in figure 22(d). Table I lists the surface 
polarization vector, the surface reciprocal lattice vector G, and the reduced surface index (which because there are 
no SSSs equals the number of zero modes per surface cell) for the seven surface cells shown in figure 14. As required, 
the count corresponds to the number obtained by direct evaluation. 

Converting Rt = 0 to Rt ^ 0 modifies the surface count by the term Rt ■ Gj{2 tt) in v. Since G has opposite 
signs on opposite parallel surfaces, the effect is to move zero surface modes between the two surfaces while keeping 
the total count fixed as required. Figure 22 plots the real part of the surface penetration wavenumber, k, of three 
representative surfaces and orientations of Rt relative to the surface normals. In figure 22(a), Rt is parallel to G of 
a vertical surface, and there are four rather than the two zero modes of the non-topological surface of figure 14 with 
surface cell 4. Thus on the opposite surface (with surface cells 5 an 6 in figure 14), there are no zero modes. Similarly 
in figure 22(b) with Rt at 30° to G of a horizontal surface, there are two zero surface modes whereas the same surface 
(cell 1 in the Rt = 0 lattice of figure 14) has only one. Finally figure 22(c) with Rt at 60° to G of a vertical surface, 
there are three rather than two surface modes per unit cell. A striking feature is that the curves versus the surface 
wavenumber q oi k for acoustic modes when Rt yf 0 is that their approach to 0 as g —>■ 0 is quadratic rather than 
linear in q. This is a consequence of the modes with q^ dispersion shown in figure 21. 


3 . Continuum limit 

As we have seen, there are two related long-wavelength properties of the topological lattices that differ from 
their non-topological counterparts: Their spectrum has peculiar low-frequency lobes [figure 21(c)], and penetration 
wavenumbers of their acoustic modes approach zero with q as q^ rather than as q. These properties must be reflected 
in the in the form of the long-wavelength elastic energy, and indeed they are. For simplicity we focus on states with 
X = {xi,X 2 ,X 2 ', 0), where a :2 > 0 is fixed and Xi is allowed to vary. The elastic energy density / can be written 

/ = ^K[{uxx - riUyyY + ^r^uly], (83) 

where ri oc xi for small Xi, while r4 > 0 and K are positive and smoothly varying near xi = 0. Thus, the Rt = 0 
and Rt yf 0 sectors are distinguished by the sign of ri. The Guest mode [49], for which / = 0, corresponds to shape 
distortions with u^x = xiUyy and Uxy = 0. When ri > 0, Uxx and Uyy have the same sign, and the distortion has a 
negative Poisson ratio [86], expanding or contracting in orthogonal directions (a feature shared by the twisted kagome 
lattice [83]); when ri < 0, u^x and Uyy have the same sign, and the distortion has the more usual positive Poisson 
ratio. Finally when ri = 0, uniaxial compressions along y costs no energy. Note that this energy consists of two 
independent positive definite quadratic forms as it must in a periodic isostatic lattice with two SSSs. Figure 24 shows 
the effects of evolution of the nonlinear Guest mode in response to compression along the cc-axis of various topological 
lattices. 

Expanding det C for small q provides useful information about the bulk- and surface-mode structure. To order 

det C = A[ql + riq^ + i-f{ql - 3qxql)] + 0{q^), (84) 

where A > 0 and 7 > 0 for small xi. The quadratic part of the equation follows from (48) and the elastic matrix, (28), 
with Kxxxx = K, Kyyyy = , and Kxyxy = 'f’ 4 K, associated with the free energy of (83). Long-wavelength zero 
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(a) (b) (c) 


FIG. 24: Evolution of the nonlinear Guest mode in topological lattices in response to uniform compression along the a;-axis. (a) 
X = {—0.05,0.05,0.05}: Small compression along x produces area change along with pure shear with a positive Poisson ratio 
in agreement with (83). At large compression, there is contraction along both x and y reflecting a nonlinear negative Poisson 
ratio, (b) X = {0.05, —0.05, 0.05}: Gompression along x now produces a simple shear component with length contraction in all 
directions (negative Poisson ratio), (c) X = {—0.1,0.025,0.075}: Gompression along x now produces some simple shear, but 
the Poisson ratio {—Uyy/u:^^) remains positive for all compressions. 


modes are solutions of det C = 0. The quadratic term, which corresponds to the elastic theory reveals an important 
difference between the bulk acoustic modes of Rt — 0 and Rt ^ 0 . In the former case, ri > 0, det C = 0 only at 
q = 0. For ri < 0, though, to order det C = 0 for = ±-\/|?’i|gy, so the elastic theory predicts lines of gapless 
bulk modes. The degeneracy is lifted by the g^ term, leading to a g^ dispersion along those lines, which can be seen 
by the cross in the density map of figure 21(c). 

As we have seen, det C(q —>■ 0) vanishes for complex wavenumbers associated with zero-frequency Rayleigh surface 
waves. Writing q — q±n -|- g||.2 x n for a surface whose outward normal n makes an angle 9 with x, there is an w = 0 
Rayleigh wave with penetration depth k~^ = |Imgj_|“^ if Img_L < 0. To order gj^ there are two solutions, 

± _ . sin6>±iyTIcos6l VG(3 + ri )7 2 /o.n 

^cos0 =F iyhq sin6**^^^ 2(cos0 ± iyhq sin0)3*^ll 

When ri > 0, the linear term is always finite and nonzero, and Im have opposite signs, indicating that there can 
be acoustic surface zero modes on all surfaces. These are the classical Rayleigh waves predicted by the elastic theory 
[59], with penetration depth 0(g|]'^). When ri < 0, the linear term in g|| is real and Hy = ImgJ oc gj^. The number 

of long wavelength surface zero modes depends on the angle of the surface. When |0| < 0^ = cot“^ ^jn], Img]^ are 
both positive, and there are no acoustic surface zero modes. The opposite surface, \6 — tt\ < 9c, has two acoustic 
surface modes. For 9c < 9 < tt — 9c, Img^ have opposite sign, so there is one mode. This is consistent with the results 
shown in figures 14 and 22. In the former figure, ri > 0, and there are acoustic zero modes on every surface. In the 
latter, Rt is chosen so that the coordinate system can always be rotated so that its direction corresponds to the —x 
direction. Thus, in figure 14(a), the surface normal G and Rt are parallel, implying by the above considerations 
that 9 = TT and that there should be two acoustic zero surface modes. In (b), 9 = Stt/G, and there are also two zero 
acoustic modes whereas in (c), 9 = and there is only one acoustic mode. This is consistent with the above 

long-wavelength analysis if tt/6 < 9c < tt/S. Finally, for surfaces such as those with cells 1, 2, 3, and 7 in figure 14, 
9 = ±7r/2 for systems with Rt in the ±x direction, and there should be one acoustic zero mode on each surface as 
there are. 


VII. REVIEW AND FUTURE DIRECTIONS 

In this review, we have presented a pedagogical introduction to Maxwell frames, in free versions of which Nb = 
dN — d{d + l)/2 and in the periodic versions of which Nb = dN, and to related frames on the verge of mechanical 
collapse. We made extensive use of the Index theorem [2], which relates the number of zero modes and SSSs of a 
frame to its site number N and bond number Nb, and of a relation between elastic energy and SSSs [47, 67] to frame 
our discussion of the elasticity and vibrational spectrum of these frames. We concentrated on Maxwell lattices, whose 
sites can be collected into unit cells whose origins lie on sites in a Bravais lattice, and we paid particular attention to 
the surface zero-modes that necessarily arise when periodic Maxwell lattices are cut to create free surfaces. 

All of the physical examples we studied were two-dimensional, both because important concepts are most easily 
explored in two rather than higher dimensions and because there is very little work of the type we discuss on higher 
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dimensional systems. There is nonetheless, interesting work to be done on three-dimensional systems. The most 
obvious frames to investigate are variants of the pyrochlore lattice, a generalization of the two-dimensional kagome 
lattice composed of corner sharing tetrahedra arranged so that lines of bonds form straight filaments. Under periodic 
boundary conditions, this lattice is a Maxwell lattice with Nb = QN. Preliminary work [82] indicates that this lattice 
can be distorted in much the same way as the kagome lattices to gap the spectrum and produce topologically distinct 
states with protected interfacial zero modes. Three-dimensional Maxwell lattices other than the simple cubic lattice 
include the many zeolite lattices, which like the pyrochlore lattice are based on corner sharing tetrahedra and which 
are in a sense 3D generalizations of the distorted kagome lattices like the twisted lattice of figure 3(d) or the more 
complex layered lattice of figure 16. These lattices all have a “flexibility window” that, like the twisted kagome lattice 
allow for easy change in volume [97]. It is likely that judicious rearrangements of lattice sites in these lattice will yield 
different topological states. 

We considered here only linearized elasticity and vibrational structure. Maxwell frames obviously have interesting 
nonlinear properties, which in the end are responsible, for example, for the ability of the kagome lattice to undergo 
large area change upon twisting. References [98, 99] explored the nonlinear properties of the one-dimensional model 
[60] discussed in section VIA and found that under appropriate boundary conditions, the surface zero-modes become 
elevated to zero-energy nonlinear topological modes that propagate freely throughout the bulk. The system is a 
nonlinear mechanical conductor whose carriers are nonlinear solitary waves not captured by the linearized theory. 
An obvious question is whether similar behavior will be found in appropriately designed two- or three-dimensional 
frames. 

As we have seen Maxwell lattices exhibit a surprisingly rich variety of vibrational responses. Ideal lattice-structures 
are generally the exception rather than the rule, and one can ask what effect do defects like dislocations have on the 
linear and nonlinear vibrational spectrum of Maxwell lattices. Reference [100] studied just this question in dislocated 
topological kagome lattices and found that zero modes can be localized at dislocations as a result of the interplay 
between between the topological dipole Rt of the lattice and the topological Burgers vector of the dislocations. Thus 
zero modes can be localized at a point by dislocations and along a a grain boundary separating two topologically 
distinct phases. Localized modes can also be created by enclosing a region of topological type B in a lattice of 
topological type A. There is certainly the potential for interesting and perhaps useful generalizations of these ideas. 
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Appendix A: States of self stress and the elastic energy 

In this appendix, we present a derivation of (33). Our starting point is (19) for the stretch energy in terms of the 
Nb dimensional bond-elongation vector E, which we break up into an affine part Eas with components given by (31) 
and a part CD = Q^U describing additional elongation that relaxes in response to the macroscopic strain imposed by 
Eaff. Multiplying on the left by the decomposition of the IVs-dimensional into the projection operators and 
P^, onto Opt. Commun. (Q) and ker(Q), respectively, and on the right by the decomposition of the dV-dimensional 
unit matrix into projections operators P(? and Pf onto Opt. Commun. (C) and ker(C), respectively, can be 
decomposed into -I- + Qsz> where = PjQ^P(?, and so on. When multiplies any vector 

to its right, the components of that vector in ker(C) are annihilated, and only the components in Opt. Commun. (C) 
survive. Thus, Q^U = -I- Q^^)Ur. When this quantity is multiplied on the left by any As-dimensional vector 

W, the term involving vanishes because it is orthogonal to Opt. Commun. (Q) and it gives zero for any vector 

in ker(Q). Thus, in Vei, Q^U can be replaced by Q^rUr- Similarly, U^Q can be replaced by U^Qrr, and the elastic 
energy as a function of U and EaS is 

U/ei = i(Efff + U^Q„)k(Eaff + Q^,U,). (Al) 

Minimizing over the dN — Nq independent components of U^. (or U^) yields 

u|^ = Q..k(Eaff + Q^,U,) (A2) 


or, because both and k^r are invertible, 


Eaff,r "b — (krr) kj-^Eaff^s 


(A3) 
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and 


Vfel = - ^sriKr) ^krs)Eaff,s = 5 Effj^[(k ^)ag] ^Eaff,s, (A4) 

where (k~^)sa is the projection of k~^ onto the null space of Q. This is (33). The final relation can be derived as 
follows: Let p the inverse of k and decompose the two matrices into their projections onto the range and nullspace of 

Q, 



Then kp = I implies that ks^pas + k^rPrs = Iss and k^sPsa + k^rPrs = 0. Thus, k^g + kgrPrsPgg^ = and k^g = 
-k„PrsPgg\ and finally kgg - kg^k'^^k^g = p”^ = [(k“i)gg]“i as required. 


Appendix B: Zeroes of det K(q) and Guest Modes. 


In this appendix we derive (47) relating the determinant of the long-wavelength dynamical matrix and the eigen¬ 
vectors of the elastic Guest mode. We begin by representing Sij as \{qiUj -\- qjUi). The the strain vector becomes 

( iqxUx \ / qx ^ \ / \ 

iqyUy = i 0 qy \ = Sa^Ui, (Bl) 

^{qxUy + qyUx)J \yy \qx) ^ 

where the summation convention is understood, a = 1,2,3, and i = x,y. The elastic energy density for a fixed q is 
thus 


fel — — o (q)rij (^), 


where = Sf^KabStj- Then using (46), 


— ^ia ( ^ ^ ^p'^pa^pb j ^bj — ^ ^ ^p^pi^pj ■; 

\p=l,2 J p=l,2 


where Spj 


Thus, 


and 


VpbSbj, or 

511 = Vi ■ {qx,0,^qy) = qxVii + ^qyVi3 

S 22 = V 2 ■ (0, qy, \qx) = qyV 22 + \qxV 2 z 

5 12 = Vi ■ (0, qy, \qx) = qyVi2 + ^qxVis 

S 2 I = V2 ■ (qx, 0, ^qy) = qxV21 + \qyV23- 

^ KiS\i+K2S2i 7^1511512-1-7^2521522^ 

^ \ 7 L 1 S 12 S 11 -I- 7 L 2 S 21 S 22 AllS 12 -I-772522 / 

detK = (77 iS?i+772S2i)(77iS?2 +7725^2)-( 771S11S12+772S21S22)" 
= 77i772(SiiS 22 — Si2S2l)^ 

= 77i772(detS)2, 


where 


det 5 = SiiS 2 j{vuV 2 j - V 2 iVij). 

The three-component vectors Vi, V 2 , and Vq form an orthonormal triad such that (vuV 2 j — V 2 iVij) 
tijk is the anti-symmetric Levi-Civita symbol, and 

det 5 = SiiS 2 je,jkVok = - ^{vo 2 ql - ^Vosqxqy + t’oWy)- 


(B2) 

(B3) 


(B4) 

(B5) 

(B6) 

(B7) 

eijkVok, where 
(B8) 


Using this equation and (B6) yields (47). 
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Appendix C: Compatibility matrix for kagome lattices 


The kagome lattice such a central role in this review. Here we derive and display the q-dependent compatibility 
matrix for a general lattice for the symmetric unit cell of figure 10 (c), which identihes and labels the three sites and 
six bonds in the cell. The labeling of the sites and bonds does not change when the cell is distorted to produce 
the twisted kagome of topological lattices. The equations relating the six bond elongations in the unit cell to sites 
displacements in that and nearest neighbor cells are 


ei = {ui{R) ~ U 3 {R)) ■ bi 
63 = (usiR) - U2iR)) ■ h 
65 = {ui{R + ai) - U 2 {R)) ■ h 


62 = {u 2 iR) - ui{R)) ■ 62 

64 = {u 3{R + ai) - ui{R)) ■ bi 

66 = (m 2(-R +ai) - M3(-R)) • ^6, (Cl) 


where R is the position of the unit cell, and ai are the lattice vectors shown in figure 10(c) and the resultant 
compatibility matrix (in a gauge in which the site and bond basis vectors and rp are set to zero) is 


bix 

bly 

0 

0 

bix 

-bly ^ 

— b2x 

-hy 

b2x 

hy 

0 

0 

0 

0 

— b3x 

-biy 

bix 

hy 

bix 

— biy 

0 

0 

e*'i‘"ib4x 


(q-^bsx 


—box 

-boy 

0 

0 

0 

0 



^Ox 

-hy ) 


(C2) 


Appendix D: Gauge choice and surface zero modes 


This appendix presents details about properties of det C in different choices of the unit cell (different gauges) and 
how switching from the symmetric unit cell (i.e., those shown in hgure 10 ) to ones that match surface boundaries 
increases the count of topological count. 

We begin with the observation that det Csym is the compatibility matrix in the symmetric gauge [(C2)] of the 
kagome lattice scan be expressed in terms of four parameters go, ...,g 3 (each a function of the unit vectors bg) as 

det Csym = 91 S 2 S 3 + g 2 *S' 3 S'i + g 3 SiS 2 + qo-S'iS' 25 ' 3 , (Dl) 

where 5'„ = w„ — 1 with w„ = Note that there is no constant term or terms linear in Wn- This is a reflection 

of the fact that the translation invariance requires the stiffness matrix K = /cC^C (43) have two eigenvalues that 
vanish as q in the long-wavelength q —>■ 0 limit and thus that det K vanish as or equivalently that det C (for any 
gauge) vanish as q^. With the help of the identity W 1 W 2 W 3 = 1, (Dl) expressed in terms of the w„s becomes 

3 3 

det Csyn,^ = G. -f ^ ( b^Wji ^ , (D2) 

n—1 n—1 


where a = gi + g 2 + 53 , K = go + g-a - a, and c„ = - go, or setting wi = z^, W 2 = Zx and W 3 = Zx ^^"^Zy \ 

where Zx = 6*^“= and Zy = 

det Csym = {a + hzx + ciz“^) -b {b 2 Z~^^'^ + 632 :^^)% + ihZx^^'^ + C 3 z]J’^)Zy'^. (D3) 

This equation can be used without further modihcation to study the topological number associated with det Csym for 
q-vectors an paths between equivalent points in the BZ parallel to the x and y axes. 

The Cauchy argument principle provides a relative count of the number of zeros and poles of any meromorphic 
function F(z) of a complex number z in the interior of a contour C: 


2 tti Jc F{z) 2111 Jc 


d\nF{z) 

dz 


dz = n — p, 


(D4) 


where n is the number of zeros and p the number of poles counted with their order (e.g., if F(z) = z p = 2) in the 
region bounded by C. If the contour C is the unit circle, the the magnitude at z at any zero of F{z) is less than unity. 
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Thus, the argument principle applied to det C with C the unit circle will count the number of zeros with \z\ < 1, and 
thus the number of decaying surface modes, minus the number of poles. 

Consider first the case for q along the y-direction (i.e. perpendicular to a surface parallel to the x-axis). Zy remains 
unchanged under the transformation qy ^ qy + An/y/Sa or equivalently under <7 —> <7 + Gi, where Gi is defined in 
figure 10 , reflecting the fact that the width of the symmetric unit cell along y is ay/3j2. det Csym is equal to Zy^ 
times a quadratic polynomial in Zy and therefore has one pole and either n = 0,1, or 2 zeros. Thus by the argument 
principle, 


1 

27ri 


d 

dZy 


In det Csyni 


= n — 1 = TTly , 


(D5) 


where rUy = —1, 0,1. In the non-topological lattices, niy = 0 indicating that as a function of qy there is one zero in 
det Csym(%) within the unit circle. In the topological lattice shown in figure 21 (c), Rt = —cti, and rUy = 1 (for the 
configuration shown in figure 21(b)) indicating two zeros in det Csym within the unit circle. Direct evaluation of the 
the integral in (D5) confirms this result. 

For <7 parallel to x [perpendicular to surface parallel to j/-axis as in figure 22(a) or (c)], det Csym is a function of 

z = and is equal to z~‘^ times a fourth order polynomial in z, and it has two poles and four zeros in z. The 
topological integral of In det Csy^ around the unit circle in z (i.e., for q^ from 0 to Itt) can, therefore, take on values of 
nix = -2,-1,..., 2. For the non-topological lattice, rux = 0 for all orientation of Rt relative to the x-axis, indicating 
two zeros within the unit circle. For the topological lattices, nix = 2 and nix = 1, respectively, for configuration figure 
22 (a) with Rt = —ai and configuration figure 22 (a) with Rt = 03 . 

Transformations from the symmetric unit cell to cells compatible with specific surfaces constitute gauge trans¬ 
formations in which the phase of detC changes by a factor of the local polarization charge of (79), —q ■ Rl = 
-dJ2rT + J2f3ri3: 


det Csur = e det Csym- (D 6 ) 

The phase factor can be determined directly from (Cl) relating bond stretches to site displacements. For example, the 
left-hand (LH) surface cell of figure 22(d) requires the displacement of site 1 through 02 and bonds 4 and 6 through 
02 and — 03 , respectively. A bond translation through a lattice vector a changes the arguments of the displacements 
of sites at the both ends of the translated bond by a. For example, a displacement of bond 4 through 02 with no 
site displacement would change the equation for the stretch of bond 4 from that in (Cl) to 64 = {u^{R -f Oi -|- 02 ) — 
Ui{R + a 2 )) ■ f> 4 -. On the other hand, a translation of a site requires the change in the arguments of the displacement 
of all occurrences of that site by the negative of the translation. Thus the equation for bond 4, which includes a 
displacement of site 1 whose argument must be augmented by the negative of that displacement, i.e., by —a 2 in the 
case of the LH surface cell of figure 22(d), becomes 64 = ( 1 * 3 ( 1 ? -f ai) — Ui{Rj) ■ 64 - as can be verified by calculating 
the expression for 64 direction from the surface cell. 

Thus for the LH cell of figure 22(d), Rl = —ai, and = z^, and detCsur = det Csym is simply 

a fourth-order polynomial in z = The two poles in det Csym have been removed, and the topological integral 

around the unit circle of z counts only the zeros of det Csym within the unit circle, four in the case of figure 22 (a) 
and three in the case of figure 21 (c). 

The right-hand (RH) surface cell of figure 22(d), on the other hand, requires only the translation of bond 6 through 
—a^, and det Csur = det Csym) which has no poles in Zy, is simply a quadratic polynomial in Zy. Again, the 
topological integral around the unit circle of Zy counts only zeros of det Csur- 

It is clear that Csur for a given surface is identical with e*'*'®'’ replaced by A, where G is the reciprocal lattice vector 
associated that surface, is identical to Csur calculate from the surface cells directly as discussed in section V C 2, (65) 
to (67). The latter equations are set up from the beginning so that there are no negative powers of A in Cgu^, so the 
integral around the unit circle of det Cgur give the count of the number of zero modes of that surface. 
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